DESY 96-048 
TUM-T31-86/96 



hep-ph/960453C] 



April 1996 



The Complete |AS| = 2-Hamiltonian in the 



Next- r 



b- Leading Order ft 



Stefan Herrlich f\ 
DESY-IfH, Platanenallee 6, D- 15738 Zeuthen, Germany 

Ulrich Nierste | 
Physik-Department, TU Miinchen, D-85747 Garching, Germany 



Abstract 

We present the complete next-to-leading order short-distance QCD cor- 
rections to the effective | AS | = 2-hamiltonian in the Standard Model. The 
calculation of the coefficient 773 is described in great detail. It involves the 
two- loop mixing of bilocal structures composed of two |AS| = 1 operators 
into |AS|=2 operators. The next-to-leading order corrections enhance 773 
by 27% to 

+0.03 

7/3 = 0.47_o.o4 

thereby affecting the phenomenology of ex sizeably. 7/3 depends on the 
physical input parameters m t , m c and A^r only weakly. The quoted error 
stems from renormalization scale dependences, which have reduced com- 
pared to the old leading log result. The known calculation of 771 and 772 is 
repeated in order to compare the structure of the three QCD coefficients. 
We further discuss some field theoretical aspects of the calculation such as 
the renormalization group equation for Green's functions with two opera- 
tor insertions and the renormalization scheme dependence caused by the 
presence of evanescent operators. 
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1 Introduction 

|AS|=2 transitions induce the mixing between the neutral Kaon states K° and 
K°. The investigation of K° — K° -mixing has revealed a lot about the short 
distance structure of nature: In 1970 Glashow, Iliopoulos and Maiani (GIM) 
postulated the existence of the charm quark |L| from the suppression of this and 
other flavour-changing neutral current (FCNC) processes. Then Gaillard and Lee 
estimated the mass of the charm quark from the measured value of the Kl — Ks - 
mass difference || . Further the violation of the CP symmetry in nature has been 
first observed in K° — K° -mixing |5| in 1964, long before the Standard Model of 
elementary particles has been constructed. The quantity ex characterizing this 
indirect CP- violation is up to now the only unambiguously determined measure of 
CP-violation in nature. Well before the discovery of the r lepton Kobayashi and 
Maskawa realized that the explanation of CP-violation within the Standard 
Model requires a third fermion family. In the subsequent decades the analysis 
of 6k has clearly been indispensable in the determination of the elements of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix. Here the CKM phase 5, which is 
the only source of CP-violation in the Standard Model, is derived as a function 
of four key parameters: the magnitudes of the CKM elements V c b and V u b, the 
non-perturbative QCD parameter B K and the top quark mass m t . ex depends on 
m t} because K° — K° -mixing is a loop process with top quarks in the intermediate 
state. As a special feature one cannot find a solution for 5 from the measured 
value of €k for too low values of the four key quantities. This has allowed to 
derive lower bounds on mt in the time before the top discovery. Yet now in the 
top era one can use the measured value for m t to constrain the allowed region for 
the CKM parameters ||. But also the accuracy of the other three parameters in 
the game has made significant progress in the last few years. To keep up with 
this progress the theorist's tools to predict the strength of the | AS| = 2 transitions 
must be sharpened as well, as we will show in the following. 
To be specific, let us look at the |AS| = 2-hamiltonian: 



16tt 



2 1V± W 



^4 + 



+ 2X c X tV3 S( 



9 2 

m c m t 



K^)Qs2{^) + h.c. 



Here Gf is the Fermi constant, Xj = VjdVL comprises the CKM-factors, and Qs2 
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is the local four quark operator (see Fig. [I]) 

Qs2 = («i7At(l - 75)^)(sfe7 M (l - 75)4) = {sd) V -A{sd)v-A (2) 

with j and k being colour indices. The Inami-Lim functions S(x) and S(x, y) || 
depend on the masses of the charm- and top quark and describe the |AS|=2 
transition amplitude in the absence of strong interactions. They are obtained by 
calculating the lowest order box diagrams depicted in Fig. ||. 




Figure 1: The diagram for the matrix 
element of Q52 to order a°. The cross 
denotes the insertion of the effective 
AS = —2 operator Qs2- 



d u,c,t s 




s u,c,t d 



Figure 2: The lowest order box dia- 
grams mediating a |AS| = 2 transition. 
The zigzag lines stand for W-bosons 
or fictitious Higgs particles. The di- 
agrams rotated by 90° must also be 
considered. 



We will be interested in the short distance QCD corrections comprised in the 
coefficients rjx, 772 and 773 with a common factor &(//) split off. They describe the 
effect of dressing the lowest order diagram in Fig. § with gluons in all possible 
ways. The 7yVs are functions of the charm and top quark masses and of the QCD 
scale parameter Aqcd- Further they depend on various renormalization scales. 
This dependence, however, is artificial, as it originates from the truncation of the 
perturbation series, and diminishes order-by-order in a s . 

The hadronic matrix element of Qs2 between the neutral Kaon states is 
parametrized as 

(K 5 I Qsiip) I K°) = \^/k<- (3) 

Here rnx and fx are mass and decay constant of the neutral K meson and \l is 
the renormalization scale at which the short distance calculation of ([!]) is matched 
with the non-perturbative evaluation of (|3]). Br in ([3D is defined in a renormal- 
ization group (RG) invariant way, because the //-dependent terms cancel when 
the physical matrix element {K° \ H^ AS<F ^ \ K°) is expressed in terms of B K . 

The first determination of (HD in the free quark model (i.e. with 77, • b(u) = 1) 
is due to Vainstein and Khriplovich and Gaillard and Lee 0. Then the QCD 
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factor 771, which is only sensitive to the first two quark families, has been calculated 
in the leading-logarithmic approximation by Vamstein, Zakharov, Novikov and 
Shifman |J . They have explicitly extracted the coefficient of the leading logarithm 
a s In (ml/M^r) from the diagrams depicted in Fig. |3] and summed this logarithm 
to all orders in perturbation theory with the help of the RG equation. In the 
same way the coefficient 772 has been obtained by Vysotskh || for the case of a 
light top quark. Then Gilman and Wise jL0) have introduced a more efficient 



method to achieve the leading log summation. Following Witten [] I] they have 
applied Wilson's operator product expansion |L2) consequently to the |AS| = 1- 
substructure and could reproduce the results of |||§ for 771 and 7? 2 . Further they 
have correctly determined 7/3, which involves a larger operator basis than 771 and 
772 due to the presence of penguin operators [[Uj. It is difficult, if not impossible, 
to achieve this calculation with the older methods of f||9|]. Further the leading 
order (LO) calculation of |TD| has only required one-loop calculations to obtain 



the leading logarithms of the diagrams in Fig. [3[ The results of |H| have later 



been extended to the case of a heavy top quark by Flynn and by Datta, Frohlich 



and Paschos [14]. 



Yet LO results suffer from certain systematic drawbacks and a precision cal- 
culation must include the next-to-leading order (NLO) terms. We sketch the 
reasons here: 

i) The fundamental QCD scale parameter Aj^g is not well-defined in the LO. 

ii) The quark mass dependence of the 7?j's is not correctly reproduced by the 
LO expressions. Especially the m t -dependent terms in 773 • S(x c , x t ) belong 
to the NLO. 

iii) Similarly the question of the definition of the quark masses (i.e. the renor- 
malization scheme and scale) to be used in (|I|) is a next-to-leading order 
issue: Hence one has to go to the NLO to know how to use m t as determined 
at Fermilab in low-energy hamiltonians like (|I]). 

iv) The LO results for 771 and 773 show a large dependence on the renormaliza- 
tion scales, at which one integrates out heavy particles. In the NLO these 
uncertainties are reduced considerably. 

v) One must go to the NLO to judge whether perturbation theory works, i.e. 
whether the radiative corrections are small. After all the corrections can 
be sizeable. 

The first step of the extension of ([I]) beyond the leading order has been done 
by Buras, Jamin and Weisz, who have derived the NLO expression for 772 [15]. 



Then we have calculated r/i in the NLO [TRl, and the present work is devoted to 



present the details of our NLO calculation of r/ 3 . For completeness we will also 
list the results of [frS],[n| for 771 and 772 and illustrate the different structure of 771, 
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F 5 



F 6 



F 7 



F 8 



Figure 3: The classes of diagrams constituting the O (a s ) -contribution to |AS| =2 
transitions in the SM; the remaining diagrams are obtained by left-right and up- 
down reflections. The curly lines denote gluons. Also QCD counterterm diagrams 
have to be included. Diagram Fq equals for zero external momenta. 



r]2 and 773. The numerical results and the phenomenological implications of our 
findings on the analysis of ex and the Kl — Kg -mass difference have already been 
given in || (for an update see [17|). We assume that the reader is familiar with 
the general concepts of Wilson's operator product expansion, the renormalization 
group and operator mixing. A detailed description of these tools in the context 
of K° — K° -mixing can be found in [|I^,|T9|j. 

The paper is organized as follows: In sect. || we first set up our notation and 
then discuss the |AS| =2 transition in the Standard Model (SM). We identify the 
large logarithms in the transition amplitude to orders a° and a] and compare 
their numerical sizes. Sect. |3| is devoted to the effective description of the | AS| = 2 
transition between the scales fi t w = 0(m t , M w ), at which the top quark and 
the W-boson are integrated out, and fi c = 0(m c ), at which the charm quark 
is removed as a dynamic degree of freedom. Here we construct the operator 
basis used in the effective lagrangian. We then match the SM amplitude to the 
effective matrix elements at the initial scale fj,tw an d determine the RG evolution 
down to n c . The latter requires the solution of RG equations for double operator 
insertions. Subsequently we describe the two-loop calculations needed to obtain 
the anomalous dimension tensor. Sect. 0] deals with the effective theory below 
the scale fi c . Here we first match the effective four-flavour theory obtained in the 
last section to an effective three-quark theory and then perform the RG evolution 
down to some hadronic scale /1. In sect. [| we summarize the analytic result and 
present an approximate formula having an accuracy of approximately 1%. Sect. |] 
contains the numerical analysis, which includes a discussion of the residual scale 
dependence of our result as well as the dependence on physical parameters. Then 
we close the paper with our conclusions. The appendices contain the results of 
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the two-loop diagrams, the renormalization factors and important RG quantities 
appearing in the calculation. 

The impatient reader may skip sects. O, 2%, 3.2.1 , 3.4.1 , 3.4.2 , 3.5.1 and 3J3 
and restrict her or his attention in sects. [| and ^| to the final results in (|122 ), 
( p5|) , O, ( P7D , Table | and the plots. 



2 |AS|=2 Transitions in the Standard Model 
2.1 Notations and Conventions 

Before writing down the result for the diagram of Fig. |^, we set up the conventions 
and notations used in this work. 

Throughout this paper we will use dimensional regularization and the MS 
renormalization scheme pD| |. Since only open fermion lines appear during the 



calculation, we can safely use a naive anticommuting 75 (NDR scheme) as jus- 
tified in |^T],[22|]. The result for 7/3 will be scheme independent, the only scheme 



dependence of ifl AS H 2 j n (TjJ) resides in the factor b(n). The scheme dependences 
of b(fi) and the hadronic matrix element in (^) must cancel, so that Bk is scheme 
independent. For the W-propagator the 't Hooft-Feynman gauge will be used, 
while the QCD gauge parameter £ is kept arbitrary. 

Let G be the | AS| =2 Standard Model Green's function, which is understood 
to be truncated, connected and Fourier-transformed into momentum space. In 
this work we are interested in the lowest order contribution to G given by the 
box diagram of Fig. |2| and the QCD radiative correction to it. The different 
contributions from the internal quarks involve different CKM factors Xj = VjdV* s , 
j = u, c, t. The GIM mechanism X t + A c + X u = allows to eliminate X u . Now the 
contributions from light internal quarks must be treated differently from those 
involving the heavy top quark. We therefore split G up as 

G = X 2 C G C + X 2 & + 2X c X t G ct . (4) 

The upper indices in the three terms in (^) denote the internal quark flavours 
involved. Further each term contains contributions with up quarks due to the 
GIM mechanism. When discussing the quark mass dependence of the G J 's we 
will frequently use the abbreviation %i = m^/M^r for the squared ratio of some 
quark mass and the W mass. In the effective hamiltonian ([I]) the three terms 
involving rji, rj2 and 7/3 emerge from G c , G l and G ct respectively. 

Frequently we will use the abbreviations L = 1 — 75 and R = 1 + 75. iV is the 
number of colours, the T a, s denote the generators of the colour group SU(N) in 
the fundamental representation, and the f abcl s are the structure constants. We 
will use 11 and U to denote colour singlet and antisinglet, i.e. (L ® R) ■ 1 means 
Sj(l — 75)4,- • Sj(l + j5)di with i, j being colour indices. The SU(N) Casimir 
factor involved will be Cf = {N 2 — 1)/(2N). We will frequently express Green's 
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functions in terms of the matrix element of the local | AS| = 2 operator Q52 defined 
in (U) and displayed in Fig. [l|. 

The &'s in (|) will be expanded in a s as 



4.71 



(5) 



The G^'^'s, n > 1, involve infrared (mass) singularities, which will be regularized 
by small quark masses m s and m^. The matrix element of some operator Q 
between quark states will be denoted by (Q) and expanded as 



(QW) = (Q) {0) 
2.2 Zeroth Order Amplitude 



+ ^M<Q( A ,)>W+0(a2( A *)). 



(6) 



In the leading order of m c /mh ea vy, where mheavy stands for m t or M^y, one can 
neglect the external momenta in (gj) and ©. 
One obtains for the three terms in (II): 



^M^S(x c ,x t )(Q S2 )M, 

lD7T z 

^- 2 M^S(x J )(Q S2 y°\ j = c,t. 



Here the Inami-Lim function || S(xj,Xk) equals 

S(xj,x k ) = S(xj, x k ) - S(xj,0) - S(x k ,0) + S(0, 0), 



(7a) 
(7b) 

(8) 



where the result of the box diagram with internal quarks j and k is denoted 
by S(xj,Xk) and the up-quark mass is set to zero. Further S(xj) = S(xj,Xj). 
Here one realizes that the effect of the GIM mechanism is not only to forbid 
FCNC's at tree level, but also to cancel the constant terms in the S"s and to 
nullify K° — K° -mixing in the case of degenerate quark masses. 

Let us look at the three contributions ( ffa|) and (|7b| ) to (f|) in more detail: 



S(x t ) 

S(x c ) 
S(x c ,x t ) 



9 1 


3 1 


3 


x t 


41 -x t 


2(l-x t )\ 


~ 2 


. 1 - x t . 



x c + 0(xl), 

—x c \nx c + x c F(x t ) + 0{x\ \vlx c 



with 



F{x t ) 



.1. I 



8xt + 4. 
4(1 -x t ) 2 



\nx t + 



x t 



4xt 



lnx 4 , (9a) 

(9b) 
(9c) 



(10) 



S 



In ( |9b"| ) and (^) we have only kept terms which are larger than those of order 
(m s m c ) I M w neglected by setting the external momenta to zero. Clearly S(x t ) is 
much larger than S(x c ) and S(x c ,x t ) reflecting the non-decoupling of the heavy 
top quark. The vanishing of S(x c ) and S(x c ,Xt) in the limit x c — > is sometimes 
called hard GIM suppression. In the imaginary part of iG, which is important 
for CP violation, the size of S(x t ) over-compensates the CKM suppression of the 
corresponding term in but the three terms are roughly of the same size. Con- 
versely the real part of iG relevant for the Kl — Kg-mass difference is dominated 
by G c and therefore insensitive to m t (see f|). 



2.3 Large Logarithms 

The Inami-Lim functions in (RI) contain logarithms of the ratios of internal particle 



masses. As we will see in sect. ^]5| the same is true for the QCD radiative 
corrections in fl5|), which in addition involve logarithms of the renormalization 
scale fi. We now discuss these logarithms in order to illustrate the effect of the 
forthcoming renormalization group (RG) improvement. 

When the product of such a logarithm with a s is large, one has to sum it to 
all orders of perturbation theory by using RG methods. Let us first investigate 
these logarithms in the zeroth order terms in @: (|9aD clearly contains no large 
logarithm because of lnxt ~ 1.5. Since (|9a|) contains no x c , the scale entering 
a s in the QCD radiative corrections to G f '^ is of the order of Mw or m t . Now 
<y s (Mw) I (47r) -lnx t ~ 1CT 2 and needs not to be summed by RG methods. We will 
come back to this point in sect. [3.6| . Yet the radiative corrections in G l contain 
the renormalization scale explicitly through hx^/Mw)- The non-perturbative 
evaluation of the matrix element in (|3|) is performed at a low hadronic scale, so 
that ln(fi/Mw) will also be large. In (pel ) we find a large logarithm | lnx c | ~ 8. 
It is 13 times larger than F(x t ). Further a s (m c )/n ■ |lnx c | = 1.0, so that the 



summation of this logarithm is indispensable. In (pb|) we would naturally expect 



the large logarithm lnx c , too. Its absence is due to the GIM mechanism, which 
we may term super-hard in this case. Yet the higher order terms in G c do contain 
lnx c , though with one power less than those in G ct . Of course G c '^ and G ct ^ n \ 
n > 1, also explicitly depend on \i. We may group the logarithms such that this 
dependence appears as ln(/i/m c ). The summation of this logarithm is performed 
by the RG evolution below the charm threshold described in sect. |j. 

The RG evolution from \x = M\y to fi = m c summing lnx c will be described 
in sect. |3|. From (|9|) one can already read off the type of summed logarithms in 
the three terms of (|j), we summarize them in Table 0. Of course there is no 
charm quark in the calculation of &, the large logarithm \nx c here emerges from 
ln/i 2 /M^/ contained in G*'^ for n > 1. If we now perform the RG evolution 
from \i = Mw down to \i — m c , we will obtain the quoted logarithm. 
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Order 


G c 


G* 


G ct 


LO 


(a s In x c ) n 


(a 3 \nx c ) n 


(a s In x c ) n ln x c 


NLO 


a s (a s lnx c ) n 


a s (a s \nx c ) n 


(a s In x c ) n 


NNLO 


a 2 (a s lnx c ) n 


a 2 (a s lnx c )™ 


a s (a s mx c ) n 


m r dependence 


none 


in LO 


in NLO 



Table 1: Logarithms summed by the RG evolution from M\y down to m c for 
the three terms in (0), n = 0,1,2,... The last line shows the order in which 
the dependence on m t enters. From the column labeled by G ct one reads off 
that the phenomenologically interesting m t dependent terms in S(x c , x t ) in ( |9c| ) 
actually belong to the NLO. This emphasizes the importance of a complete NLO 
calculation for t] 3 . 



2.4 The Definition of Quark Masses 

When discussing analytical expressions beyond the LO, one must specify the 
definition of the quark masses. This point is often handled incorrectly in phe- 
nomenological analyses, so that we discuss it in some detail now. 

Any perturbatively calculated interacting fermion propagator is proportional 

to 

i. 

(ID 



p — m + £ (p 2 , m) 

Here m is the renormalized current fermion mass, which enters the Lagrangian, 
and iS(p 2 ,m) is the 1PI self-energy describing the dressing of the free fermion 
propagator. £ starts at second order in the gauge coupling g and may be cal- 
culated to some order g 2n . Now different renormalization schemes may involve 
definitions m and ml of the fermion mass, which differ by a perturbative series: 



m = m 



(l + Kg 2 + 



Yet also E in (11) is different in both schemes, but the position of the pole in 



LTD is the same within the calculated order: 

' -£' (m' 2 ,m') = m-S(m 2 ,m) + O (g 2n+2 ) 



m 



The freedom in the choice of the mass counterterms allows us to move any desired 
constant term from m to S. If the fermion is a lepton and therefore exists as 
a free particle, m is commonly defined as the pole mass m po i e corresponding to 



m poie' m poie) = 0. Since the pole at p 2 = m 2 le in (jTTJ) is observable for free 



fermions, m po i e is sometimes called the physical mass. Yet the strong interaction 
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confines quarks into hadrons and the quark pole mass is not observable. In fact 
the infrared structure of QCD imposes a strongly divergent perturbation series 
upon observables expressed in terms of m po i e , which is most likely only a suitable 
parameter for very low orders of perturbation theory [p3 ]. Instead in QCD one 
preferably uses a short distance mass such as the running quark mass m(/x) in 
the MS scheme. It has the additional advantage to allow for a simpler solution 
of the RG equations. Its relation to the one-loop pole mass reads: 



m 



(i) 

pole 



m(/x) 



l + ^C>f4 + 31n4 N 



47T 



(12) 



Clearly the proper definition of the quark mass only matters beyond the leading 
order. The one-loop relation ([12]) is the appropriate one for the NLO calculation 
presented in this paper. At Fermilab the top quark pole mass is measured. w^p 0le 
is larger than mt(mt) by a factor of 1.045 corresponding to 7-8 GeV. 

Finally, if the renormalization scale /i is much different from m, one must 
sum the logarithm in ( |T2] ) to all orders with the help of the RG, see 
appendix |C[ 



in 



2.5 The 0(a.) Corrections 



The 0(a s ) corrections to the box diagram (see Fig. §) were first evaluated in (15 



for the case of arbitrary internal quark masses. These corrections have been 



necessary to obtain r/i and r/2 in the NLO (T^,|T6[. We stress here that one does 
not need them for the NLO calculation of r/3, which is the novel issue presented in 
this work. Nevertheless it is instructive to look at G ct ' W as well for three reasons: 
First one can identify the logarithms summed by the RG evolution, which provides 
a very good check of the results presented in sects. |3| and f|. Second one can partly 
estimate the size of the next-to-next-to-leading order (NNLO) terms. Third the 
0(a s ) terms will be useful in the discussion of the proper treatment of the physics 
between the scales \i = m t and /1 = Mw presented in sect. [3~6 . 
Generally the 0(a s ) terms are of the form 



l07T z 1 



Here new operators have emerged^] 



>( ) + /4(T>(°) + /^> (0) } 



T 



U 



(L 
1 



L + R® R - a 



[IV 



a 



fj,v • 



( 7m L ® ^R + lt ,R ® YL) 



^±1 

2N 
'N 2 + N 



j = c,t,ct. (13) 



(14a) 



2N 



-1 - 



2 AT 



(L&R + R&L) 



N 2 + N - 1 
2iV 



11 - 



2N~ 



(14b) 



^We omit the spinors on the external quark lines. 
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which are written in a manifestly Fierz self-conjugate way. In the following we 



will discuss the coefficient functions h in fll3|) in great detail, starting with those 
of the new operators: 



h T S(x c ,x t ) 



h v = huS(x c ,x t ) 



hn 



-3-C; 



h 3 T 



h. 



h T S(xj), 
huS(xj), 



J 
J 



c,t; 
c,t; 



3 + £ m d m s m\ 



mi 



mt 



(15a) 
(15b) 

(15c) 



We first observe that flTcf ) is obviously unphysical, because the functions in QT5D are 
gauge dependent. This is an artifact of the use of small quark masses to regularize 
the infrared singularities while at the same time using on-shell quarks with zero 
four-momentum for the external states. For the same reason we encounter the 
new operators T and U. Yet the one-loop matrix element of Qs2 corresponding 
to the Feynman diagrams of Fig. [| involves the same operators (14) with the 






Figure 4: The classes of diagrams constituting the O {a ^-contribution to the 
matrix element of Qs2', the remaining diagrams are obtained by mirror reflec- 
tions. The cross denotes the insertion of Qs2, the curly lines gluons. A QCD 
counterterm has to be included. 



same coefficients ||T5|,|l6t|21 



(&a(A*)> = (Qs2^)) 



(o) 



Atx 



a(/i)(Q S2 ) {0) + h T (f)^ + h v (U)W] . (16) 



We will discuss the third coefficient a(/i) in conjunction with h? in (0). Now 
(|i~3l) and (jl6|) allow to express iG^ as 



iG ct (n) 



(~<2 



167T 2 W 
(~<2 



S{x^)) + ^V{ Xj ^)^) 



16tt 2 w 
+0(« s 2 ), 



47T 



S (x c (n),x t (ii)) + 



(QS2(»)) 

j = c,t, (17a) 

(Ga»(A*)> 
(17b) 



^-k ct (x c (^x t ^)^) 
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where the new coefficient k? is related to h? in flT3| ) and a(/i) in (O) via 

h? = k 3 (xj, fi) + a(fi)S(xj), j = c,t, (18a) 
h ct = k ct (x c , x t , fi) + a(n)S(x c , x t ). (18b) 



Now the unphysical terms hx and hu containing £ and the infrared regulators 
nid and m s have been absorbed into (Qs , 2) ( ' 1 ' ) in (0). Likewise these unphysical 
terms in h? have gone into a(/z) in (|T8|), so that A;- 7 only depends on x c , x t and /x. 
In the coefficient a(/z) in ( |l6l) we also split off the gauge and IR parts: 



N — 1 mdm s 
a(fi) = cH „ , T 6m- 



+2£ 
-3 C F - 5 



2iV 



N — 1 



2N 

N - 1 
2iV ' 



1 - 



ra 2 In -f - m 2 . In -/ \ jV — 1 _ m d m 



m s — TO 2 , 



+ 



2iV 



In 



(19) 



Next we write down the results for the k J, s grouped according to powers of the 
large logarithm lnx c : 



k c (x c ,x t ,fi) 



-3x c In 2 x c + x c In x c I -5 - 6 In — + kf ing a (x t ) 
t V 2 ct 

JrX cKing ( X t) + x c hi 2^sing,a ( X t) 
TTXi n 



N-l 

2N 



-9x c In x c + x c \nx c [3 + k c oct a (x t ) - 6 In 



+ x cK t ct (x t ) + x c lia —kZt (x t 
mi 



+c x c In x c — c x c F(x t ) + O {x 2 c In 2 x^j 



mt 



(20a) 



k°(x c ,{i) 



k\x t ,fi) 



—x c + 6x c In 

mi 



N - 1 
2N 



u 2 4vr 2 

12x c In x c + 6z c In — - — llx c H x c 

to 2 . 3 



fcL„ (z t ) + 6 In -^-fe* 



stng 



N-l 
2N 



i\/r2 '"smg,< 
1V1 W 



kl ct (x t )+6\n-^S(x t ) 
lvl w 



cS(x t ). 



(20b) 



(20c) 
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We have hidden a complicated x t dependence in the following functions: 

K s ing{X t ) ~ 6X t ^ _ 2 + 7T ^ 

2t . /i \ -2-2xt + x 2 



+3x t Li 2 (1 - x t ) — — -3— 

1 - x t ) 

20 - 12x t - 51x, 2 - llx 3 + 6xf 



2-2x t + x 2 

2{\-x t y 



+ lnx t - 



4(i-^r 



2 4 — 12x t — x\ — x\ — 3xj + 



+ 31n<s t (21a) 



., , , _ - 32 + 8x t + 3x 2 | ^ 2 -8-7x 2 + x 3 
4x t (x t - 1) 6x 2 



*=<<*> = 4s.(,.-l> + " 
+ Li 2 (l-x t ) 



8 - 16x t + 23x 2 - 29x 3 + 4x? + x t 5 



(x t - ly 

-39 4- 3fi<r. - 1 9^ 2 - 
+ lnx t - 



-32 + 36x t - 12x 2 t - 7x 3 



4x t (x t - 1)' 



2 12 - 16x t + 5x 2 + 2x 3 

+ ln x t — -2 , (21b) 

4(x t - 1) 

, t / \ 4 - 39x t + 168x. 2 + llx 3 n n T . . 5 + x t 

= *t 77 773 ^ + 3x 3 Li 2 (l-x t )- ' 

4 (xt - 1) {x t - 1) 

-4 + 24x t - 36x 2 - 7x 3 - x\ 

+3x t lnx t — -4 

2 (x t - 1) 

, o 3 1 2 13 + 4x t + x 2 

+3x; i ln x t — - —T-, (21c) 

2 ^ - 1) 



-64 + 68x t + 17x 2 -llx 3 , 8 
" 4(1^ + 3^ 



N 8 - 24x t + 20x 2 - x, 3 + 7xf - x, 5 
+2Li 2 (l-x t ) ' 



x t (x t - 1) 
-32 + 68x t - 32x 2 + 28x 3 - 3x^ 
+ ln *> 2 (*.- !)' 

+ ^, ^|i±MzM , (21d) 

2 (x t - 1) 
x t , -4 + 12x t - 3x, 2 + x? 



= -97—^ + 3 In xt 



(xt-iy 2( Xt -iy 

6fF(xt)-l + Xt^-F(xt)V (21e) 
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„2 



CW = I^T + S 4 ,^^ lnx f = 6F(x t ) (21f) 
2 — 1 2 (1 — x t ) 



, t -4 + 18a; t + 3x 2 + 

WW - 4(x 4 -l) 3 2 (x t - l) 4 



Here Li2 (x) denotes the dilogarithm function 

Li 2 (x) = -/ X ^ ln(1 - Xt) (22) 
Jo t 

and F(x t ) has been defined in fllPl). Let us now look at the ingredients of fllTD in 
more detail: ( |I7[ ) is an operator product expansion (OPE) of the Standard Model 
amplitude in terms of the local |AS|=2 operator Qs2- The terms in brackets 
are the corresponding Wilson coefficients, yet in ordinary perturbation theory 
without any RG improvement. From (Q) and (|20|) one verifies that they are 
gauge-independent and free of the infrared regulators m s and ma- If we had used 
the dimensional method also to regularize the IR singularities, the operators T 



and U and the gauge dependence would be absent on both sides of (|17D , but the 
fc J 's would be unchanged. Further the Wilson coefficients do not depend on the 
choice of the external states used in the calculation of the matrix element. Now 
fc*(xt,/x) is simply the 0(a s ) part of the initial condition for the RG improved 
Wilson coefficient needed for the calculation of r] 2 . &'(x 4 , //) is called D(x t ) in fl"5| . 
The RG evolution from /i = M\y down to a low hadronic scale sums a s ln(/x/M^) 
to all orders in perturbation theory. The situation would be the same with k ct 
and k c in a fictitious world in which the charm quark is so heavy that lnx c is 
small. To describe the real nature, however, we must first sum a s lnx c to all 
orders as well. Since this is the purpose of the subsequent sections, we discuss 
the powers of In x c term-by-term now. Therefore we have arranged ( pD| ) such that 
large logarithms can easily be distinguished from small terms. 

k ct : In ( |20a|) one immediately observes two terms oc In 2 x c , which we have ex- 
pected from the fact that S(x c , x t ) in (pcf) already contains the logarithm 
lnx c . They all belong to the LO of RG improved perturbation theory, 
c.f. Table ||. Further k ct exhibits In 1 x c terms. They are linked to the In x c 
term of S(x c , Xt) and constitute the first terms of the NLO expression. Note 
that the LO terms, the one cx lnx c in S(x c ,x t ) and the ln 2 x c term of k ct , 
are independent of m t . Top dependence first enters through the NLO terms 
In x c in (Pc| ) and the functions kf ing a and kf ct a . Finally k ct contains a In x c 
piece, which already belongs to the next-to-next-to- leading order (NNLO). 
Therefore we will not need it in our analysis. 

k l \ (|20c|) contains the logarithm ln/i 2 /M^, which gets large when \x gets small. 
We expect the appearance of just this logarithm, because there are no other 
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largely separated mass scales in G 1 . The whole S(xt) in ([5a|) and the term 
oc In /i 2 /M^ in belong to the LO. The In /i 2 /M"w part of (gOc^) is a 

NLO term. 

fc c : In ( [20 b| ) we find the large logarithm lnx c , which together with S(x c ) belongs 
to the LO terms in G c . Note that the bracket proportional to Cf does not 
contain such a logarithm, although the analogous term of k l contains one. 
This term is connected with the running of the corresponding quark mass. 
This mass is small in G c but large in G l . The non-logarithmic part of ( |20 B ) 
again belongs to the NLO. 

To gain an impression of the relevance of a calculation beyond the LO, we 
further look at the numerical sizes of the S- and fc-functions. For typical values of 
the input parameters we obtain the numbers summarized in Table |2|. Generally 











a s (m c ) y 




j 


In 1 


In 


In 2 


In 1 


In 


ct 


0.002176 


0.000311 


-0.002766 


0.000203 


-0.000235 


c 





0.000264 





-0.000192 


0.000029 


t 





6.18 





-9.21 


-3.16 



Table 2: Numerical values of different contributions to the full SM amplitude. The 
listed values correspond to the case that ot s and m t are evolved down to /i = m c 
without integrating out heavy degrees of freedom. We have split the S- and k- 
functions according to the power of In x c involved. We use a[^ =6 ^(m c ) = 0.277 
and m^ -6 ^ (m c ) = 309 GeV corresponding to f — 6 active flavours and typical 
values for a s (Mz) and m t (m t ). The a] ln n x c -parts are of the same size as the 
ot° s ln ra_1 x c -terms illustrating the need to sum a s In x c to all orders in perturbation 
theory as described in sect. ||. 



the a] ln n+1 , n — 0, 1, terms are linked to the a° ln n terms via the RG equation, 
cf. Table |l|. In Table |2] one observes that the size of the a] ln n+1 contribution is 
about as large as the a° In™ one. This emphasizes the need for the summation of 
large logarithms to all orders by means of an operator product expansion (OPE) 
and RG techniques. The thereby improved result will contain the coefficients 
of the logarithms listed in Table || evaluated for ji = jj, t w = 0(m t , Mw)- For 
example one finds for fi = m t (m t ) = 167 GeV: 

— S(x c ,x t ) = -lnx c + 0.59 (23a) 

x c 
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s -^k ct (x c ,x t ,a) = 0.026 In 2 x c - 0.22 In x c - 0.27. (23b) 

x c A-k 

The large magnitudes of the NLO coefficients 0.59 and —0.22 compared to the 
LO terms further emphasize the importance of the NLO calculation. Finally 
the constant term —0.27 enters the initial condition of the NNLO calculation. 
It amounts to 46% of the corresponding NLO term 0.59. The discussed initial 
condition, however, has a much smaller impact on the complete NLO result for 
773 than the operator mixing worked out in the following section. 



3 Effective |AS| = 2 Transitions above the Charm 
Threshold 

In this section we will sum the large logarithm In x c found in ( |9"c| ) and ( ^HQ to all 
orders in perturbation theory. This is done in two steps: First one sets up an 
effective lagrangian £ e fj in which the W-boson and the top quark are removed 
as dynamic degrees of freedom. In £ e g the |AS| = 1 and |AS|=2 transitions are 
described by local four-quark operators, which are multiplied by Wilson coeffi- 
cients. The logarithm lnx c is thereby split as lnx c = ln(/i 2 /M^) + ln(m 2 //z 2 ). 
Here the former term resides in the Wilson coefficients, which are functions of m t , 
Mw and /i, and the latter is contained in the matrix elements of the four-quark 
operators depending only on /i and the light mass parameters. The second step is 
the application of the RG to the Wilson coefficients. For /i = fi tw = 0(M W , m t ) 
there is no large logarithm in the Wilson coefficients. The RG evolution from /j, t w 
down to fi c = 0(m c ) sums ln(/j 2 //j 2 jy) to all orders in perturbation theory. The 
RG improved coefficients finally multiply matrix elements which do not contain 
large logarithms, because ln(m 2 //z 2 ) is small. 

When passing with \i below /i c we must also integrate out the c-quark field. 
This will be described in sect. |j. 



3.1 General Structure of the Effective Lagrangian 

After integrating out the top quark and the W-boson we are left with an effective 
five-flavour theory described by a lagrangian of the generic form 

^tff^ 2 = 7=^CKM ^ CfcQfc ~^CKM ^ CjQj. (24) 

V2 k Z l 

Here the Vckm denotes products of CKM elements. The Qk, Qi represent lo- 
cal |AS| = 1 and |AS| = 2 operators and the Ck, C\ are the corresponding Wilson 
coefficient functions with Fermi's constant factored out. The |AS| = 1 operators 
Q k are necessary for the proper treatment of |AS| = 2 transitions, because they 
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contribute to the transition amplitude through Green's functions with two oper- 
ator insertions. An example is shown in Fig. [5|, which is simply obtained from 
Fig. |2| by shrinking the W-boson lines to a point. The |AS|=2 operators Qk 




Figure 5: Diagram Dq in the effective five- and four-quark theory. The crosses 
denote insertions of local |AS| = 1 current- current operators. 



(see e.g. Fig. [1]) can likewise be obtained by shrinking the whole box function 
with internal top quarks in Fig. || to a point. Yet the Q^s are also needed for the 
light quark contributions. Diagrams of the type in Fig. are in general divergent 
and require counterterms (omitted in (p4[)) proportional to |AS| = 2 operators. In 
general both |AS| = 1 and |AS|=2 operators in fl2~4]) contribute to |AS| = 2 tran- 
sitions. Yet there may be special cases in which either the former or the latter 
are absent. As we will see later, all three possibilities are realized in rji, r]2 and 

In the following sections the detailed structure of C^^ 2 will be worked out. 
This requires the following steps: 



Find the minimal operator basis to be used in (24) sufficient to describe 
the physics of the |AS|=2 transition. Here one must first find a set of 
operators closing under renormalization. Subsequently one can eliminate a 
set of unphysical operators. 



ii) Match the full SM Green's function G of (f§), (§), (0) and ([17]) to the one 



obtained in the effective theory and thereby determine the Wilson coefficient 
functions Ck and Cj at the initial scale fi = n t w = 0(Mw,mt)- 

iii) Next prepare for the RG evolution of the Wilson coefficients from /x = n t w 
down to the final scale // = \i c = 0(m c ). For this one must derive the general 
RG equation for Green's functions with double insertions (see Fig. ||) and 
its solution. The RG equation involves an anomalous dimension tensor in 
addition to the familiar anomalous dimension matrices. 

iv) Determine the anomalous dimension tensor in the NLO for the operator 
basis at hand. This requires the calculation of two-loop diagrams. 

Finally we discuss the size of the remaining non-summed logarithm In m t /My\r. 
We do not sum this logarithm because we simultaneously integrate out the top 
quark and the W-boson. 
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3.2 The Operator Basis 



At first we restrict the set of operators in ([24]) to the lowest contributing di- 
mension, which means dimension six for the Qfc's. As for the Qi's, we must 
distinguish whether they correspond to the SM graphs with internal top quarks 
or whether they enter C eS as counterterms in the light quark sector. In the 
former case there is only one physical |AS| = 2 operator Qs2, introduced in (|2|), 
with a dimension-two Wilson coefficient Cgl containing all information on m t and 
My/ 1 15[| . The latter operators have the same dimension as the diagrams they 



renormalize, which is eight as can be easily seen from Fig. [5] Higher dimension 
operators correspond to terms suppressed by powers of m\ ght /m\ eiwy) which we 
already neglected in the SM amplitude, see @ and (f20|). 

We will now establish the |AS| = 1 part of the operator basis, i.e. the Q^'s in 
(|24|). Consider first the SM |AS| = 1 transition of Fig. || with only light quarks 
k, I = u, c on the external legs. Contracting the W-boson propagator to a point 
yields the diagram of Fig. [7], in which the cross denotes the insertion of the | AS| = 1 
current-current operator 



Q2 1 



(s lfl Lk) ■ (j-fLd) ■ 1, 



(25) 



which we have already met in Fig. E. It is well-known that QCD corrections to 





Figure 6: A |AS| 
standard model. 



1 transition in the 



Figure 7: The matrix element of a 
| AS | = 1 current-current operator Qf 2 
displayed as the cross. 



Q2 induce counterterms proportional to other operators, so that Q 2 mixes with 
them under renormalization. 

In the case of A; ^ I the mixing is particularly simple, only mixes with 

Q\ l = (s^Lk) ■ (l^Ld) ■ 1. (26) 

The one-loop mixing proceeds through the diagrams of Fig. || Hence the 
corresponding part of C^ S ^ 2 reads 

E £ V: s V ld C^Qf h -. (27) 

V ^ i,j=l k,l=u,c 
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Figure 8: The classes of diagrams yielding the matrix element of the |AS| = 1 
current-current operators Qf 2 as we ^ as the mixing of the current-current oper- 
ators among themselves to order a s . 



Here and in the following the superscript "bare" denotes unrenormalized opera- 
tors, while renormalized ones do not carry an additional superscript. The 2x2 
renormalization matrix is diagonal in the basis 

Q± = \(Q¥±Qt l ), (28) 



provided one preserves Fierz symmetry in the renormalization process f2T|,|33 

For k = I additional operators enter the scene, the so-called |AS| = 1 penguin 
operators, which appear in different species. In [^4],|25| the NLO mixing of Qi,2 
with quark-foot penguin operators Q3 to Qq displayed in Fig. |9| has been worked 
out. They read 



Qs = 




E 

q=d,u,s,. 


(q^Lq) 


1 


Q4 = 




E 

q=d,u,s,. 


(q^Lq) 


2 


Q 5 = 




E 

q=d,u,s,. 


(q^Rq) 


1 


Q 6 = 


(sj^Ld) 


E 


(qYRq) 


5 



(29) 

q=d,u,s,... 

and enter fl24|) with Vckm = — A*. The summation runs over all active flavours, 
at present q = d,u, s,c and b. 

Yet we may doubt, whether Qi — Q 6 are sufficient to describe the |AS| = 1 
substructure in |AS|=2 transition. Indeed, Qi — Q e mix via the diagrams of 
Fig. |TT] into operators containing only two quark lines such as the gluon-foot 
penguin operators Q 9 i, Q 9 2 and Q 9 3 depicted in Fig. [12]. Likewise loop diagrams 
with Q g i — Q g 3 require counterterms proportional to Q3 — Qq (see Fig. [lCj) and 
similarly to a ghost-foot penguin operator 

Q FP = S-i^d ■ (d^rf) n c f abc . (30) 

Here r\ c denotes a Faddeev-Popov ghost field. We can easily construct a |AS| =2 
diagram with one of these operators and Q2, see e.g. Fig. [T3|. 
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Figure 9: The four-quark penguin op- Figure 10: Diagrams for the mixing of 
erators Q 3 ,.. .,6 Q 9 i denoted by the shaded cross into 

four-quark penguin operators. 



d s d s d s 




Figure 11: The mixing of a |AS| = 1 current-current operator, denoted by the 
cross, into gluon-penguin operators. Diagrams with the gluon lines permuted 
compared to G2 and G3 must be included, too. 




Figure 12: The |AS| = 1 gluon-foot penguin operators 




Figure 13: |AS| = 2 diagram with one insertion of Qi,2 ; denoted by the light cross, 
and Q 9 2, denoted by the shaded cross, which could contribute to 773. 
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Fortunately these extra operators Q g i — Q g 3, Qfp combine to unphysical oper- 
ators, which can be dropped from the renormalized effective lagrangian |26|-p0|]. 
Hence the |AS| = 1 operators in (pip can be restricted to Qf , . . . ,Qq and the 
|AS| = 1 part of £gg S ^ reads 



£|AS|=1 



G f 



V2 



E E vzv u p i ctf-\ t Y t c i Q i 

i=l k,l=u,c i=3 



(31) 



We will illustrate the irrelevance of Q g i — Q g 3, Qfp in the following section. This, 
however, first requires the determination of the local |AS|=2 operators Qi in 

(HD- 

Consider first the simplest case, which is realized in &: the SM result for & 
is proportional to the square of a heavy mass m t , Myy times a function of their 
ratios (see ([TBI), (|I3|) with j = t). A diagram in the effective five-flavour theory 
containing two insertions of | AS| = 1 operators can only produce a result involving 
light mass parameters, which we already neglected in ([7T5p, (|T3"D. Therefore the 
diagrams with double insertions do not contribute in this case and we are left 
with a |AS| =2 contribution to C^ S]j=2 of the form 

G 2 

^KG S2 Qs2- (32) 



C S 2 has been obtained in the NLO in |15 



In the light quark sector the situation is a bit more involved: Consider first the 
diagram of Fig. [5] with two internal charm quarks and zero external momenta. The 
only physical dimension-eight |AS| =2 operators required to absorb its divergence 
reads 

2 2 
— 777 ~ 777 

which follows from power counting and the absence of any non-zero mass param- 
eter apart from m c . The inverse powers of g are introduced for later convenience 
as in |T0 . One may arbitrarily shift such factors from the Wilson coefficient into 



the definition of the operator. The factor \i 2e stems from gbarc = Z g gfi e and the 
fact that 



2 

777 

gbare = _J|baxe . ^^axe _ ^ 

fi'bare 



must be independent of fJ^\. Any other dimension-eight |AS| = 2 operator contains 
one or two powers of m c less than Qi and derivatives and/or gluon fields instead. 

2 ^Here and in the following the quark fields like s and d in ( |33| ) and ( p^ ) are bare fields. 
The wave function renormalization constant is taken into account when calculating matrix 
elements 
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Their on-shell matrix elements are suppressed by powers of m s /m c with respect 
to those of Q 7 , so that they do not contribute to the coefficient of the leading 
dimension-six operator below the charm threshold (cf. ([1])). Likewise they cannot 
mix with Q7 under renormalization. 

Next we determine the CKM factor multiplying C7Q7: Generally we could 
expect a term proportional to A^ and one proportional to X c X t corresponding to 
G c and G ct respectively. Due to the special structure of G c , which does not contain 
a large logarithm to order a° in the SM amplitude (|9l5D , no term proportional to 
occurs. 

Hence J^^f^ 2 is found as 



|AS|=2 

off 



G p 



r" 1 

\2ri{t) 7-1 Abare 

Gl 



E^ 1 E Vk S VidQf' haTC - A t E z^Qf™ 

j=l k,l=u,c j=3 



-A c At 



E E CkQZkij + CV-Z77 1 Q7 

Lfe=l Z=l 

+counterterms proportional to unphysical operators. (35) 



bare 



Here the first line is the pure |AS| = 1 part (RTI) written in terms of bare oper- 
ators and the corresponding renormalization factors. The second line contains 
everything related to a single insertion of the local |AS| = 2 operator Qs2- The 
third line consists of two parts: the counterterm involving ZZ 7 renormalizes the 



matrix elements of the type in Fig. [| or Fig. [14] with one insertion of Qi and one 
of Qk- The part of ( |35D involving Z^j 1 expresses the renormalization of the local 
I AS I =2 operator Q7 in single insertions to cancel the divergences of the diagrams 



depicted in Fig. [|. The absence of a local 
3BT) is due to the GIM mechanism 



I AS I = 2 operator proportional to \ c in 



sec 



10, 161). The structure of U35D becomes 



much more transparent, if one regards the Z-f actors as renormalization factors of 
the effective couplings Ci, Cj rather than of the operators. The bracket in the 
third line of ( |3~5"1) may then simply be interpreted as the "bare" Wilson coefficient 
of Q 7 in analogy to the bare coupling in QCD. It is renormalized both by QCD 
and by effective | AS| = 1 interactions. The last line of (|35|) contains counterterms 
proportional to unphysical operators, which are described in the next section. 
One automatically includes them in the calculation, if one simply subtracts the 
sub divergences diagram-by- diagram in the effective theory. 



3.2.1 Unphysical Operators 

We will discuss three types of unphysical operators: BRS- exact operators, those 
which vanish by a field equation of motion (EOM) and evanescent operators. 
Therefore let class B contain all BRS-exact operators, class 8 the operators van- 
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ishing by a field equation of motion and class V the physical and evanescent 
operators. 

Let us first reduce the operator basis with respect to classes B and £. The 
techniques to do this have been worked out in p6|-|3l|]. They are widely used for 
the treatment of Green's functions with single operator insertions. Here the new 
issue is the application of the theorems concerning double insertions |28|,[5(J in a 
concrete calculation. 



The mixing of the classes introduced above follows the general pattern 



z- 1 



^ * * * \ 

* * 
v * j 



corresponding to class 



B 



(36) 



The block-triangular form of fl5o] ) ensures that the Wilson coefficients of the 
operators from classes B and £ do not mix into the ones of physical and evanescent 
operators in V. 

Further we know that the on-shell matrix elements of operators in class B or 
£ vanish. This is important for the matching of transition amplitudes in different 
theories, e.g. of the full SM and of an effective five-flavour theory. The vanishing 
of the on-shell matrix elements of operators in B and £ ensures that they do not 
contribute to physical Wilson coefficients at the matching scale. 

Since operators from B and £ contribute neither to the matching nor to the 
mixing of class V, we may neglect them in the discussion of the RG equation and 
evolution. 

Let us now organize the penguin zoo according to the classes V, B and £. 



Due to the theorems of Pq -|31f the operators Q g \ — Q g s, Qfp can be arranged to 
appear in the combinations 



and 



with 



Q 



gEOM 



1 



Qbrs — Qfp + Qgf £ B 



Qg~\ (<?4 + Qe) + j^r (Qs + Qb) - Qbrs 



e £ 



Q g = -s lfl LT a d-D a u F^, Q g 
9 



--s 1 .LT a d-d tl d u A a v . 



(37) 



(38) 



(39) 



The latter is stemming from the gauge fixing part of the QCD Lagrangian. Q g EOM 
vanishes by the equation of motion of the gluon field. Qbrs and Q g EOM are 
discussed in detail in [jSOf] . 

Since now Q g i — Q g3 , Qfp have been traded for linear combination belonging 
to classes £ and B, one may drop them from the renormalized operator basis 
when calculating |AS| = 1 transitions. 
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In the case of double insertions the situation is more complex, because now the 
operators of class S may give a nonzero contribution in on-shell matrix elements 
and therefore their presence cannot be ignored for the matching. Yet it is possible 
to absorb the effects of these operators into the coefficient of a |AS| = 2 operator 



30|, 31|. Such non- vanishing matrix elements with two insertions of Q g EOM appear 
in G 1 of (H). Since the effective five- or four-flavour theory does not contain m t 
and M w anymore, this contribution is suppressed by a factor of ^n g ht/ m heavy 
compared to the term in the second line of (^) and can therefore be neglected. 
Yet in the calculation of G ct we will face the operator 

QqEOM = m 2 c sR0d, (40) 

which vanishes by the quark equation of motion. Its effect on | AS| = 2 amplitudes 
can likewise be absorbed into C 7 . 

On-shell matrix elements involving one or two operators from class B still 
vanish, we can therefore drop them in the case of double insertions, too. 

Another important class of unphysical operators are the evanescent operators 
also contained in V, which generally appear in theories with four-fermion inter- 
actions, if one uses dimensional regularization. As an example one may look at 
the 0(a s ) corrections to the matrix element of Q\ t 2, displayed in Fig. || When 
calculating diagram C 2 , C 3 one faces the structure 

Q' = l P l a l^l-ls)®l tl l a l p (l-ls) = ^-ae)Q + E 1 iQ}+0(e 2 ) 1 (41) 

where Q = ^^L ® ^L. is evanescent, i.e. it vanishes for D — 4. a is an 

arbitrary real parameter, its choice belongs to the definition (|41~|) of 

When perturbative results are improved by means of the OPE and RG tech- 
niques, subtleties arise: Evanescent operators can affect the matching proce- 
dure and the operator mixing [[3l],|32|]. In J2j], 31] a finite renormalization of 



the evanescent operators has been proposed to render their matrix elements zero. 
Doing so the Wilson coefficients of the evanescent operators become irrelevant 
at the matching scale. If this should hold at any other scale, one has to ensure 
that the Wilson coefficients of the evanescent operators do not mix into the ones 
of the physical operators. This has been proven in jEJ for a very special and 
calculationally inconvenient definition of the evanescent operators. In |33j we 
achieved the following improvements: 



i) We have generalized the proof of [32] to an arbitrary definition of the evanes- 



cent operators, which includes the one chosen in |2l] 



ii) We have shown that the arbitrariness in the definition of the evanescent 
operators displayed in ( [4~T|) introduces a scheme dependence into the phys- 
ical Wilson coefficients at the matching scale as well as into the physical 
anomalous dimension matrix starting in the NLO. This distinguishes the 
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evanescent operators from the operators in classes 8 and B. Of course this 
scheme dependence cancels in the product of Wilson coefficients and ma- 
trix elements at any scale. We give explicit formulae to transform Wilson 
coefficients or anomalous dimension matrices from one scheme to another. 
These formulae are particularly necessary if one wants to combine Wilson 
coefficients and anomalous dimension matrices calculated with different def- 
initions of the evanescent operators. 

iii) We have extended the findings from the case of single insertions to double 
insertions, which is needed for this work. 

It is important to note that the first and third point above enables us to use the 



results of 21, 24, 25 for the IASI = 1 hamiltonian in 



The physical operators in V needed in our calculation are 

Qf, i = l,2, k,l = u,c; Q j: j = 3,... 6; Q 7 . (42) 

For the evanescent operators appearing in the calculation we use the definition: 



Ei[Q 3 ] = te^T'TYi-^ + a^T.^'fPij, j = 1, ... 4, (43a) 
Ei[Qj] = [wvlr 1 R ®7 r 'Y^L-(l6 + a 2 e)^R®^L}K lj , j = 5,6, (43b) 



Ei[Q 7 ] 
E 2 [Q 7 ] 



9 2 

2 

mi 



[l^lulrjL ® 7 V7 M £ - (4 + 5ie) j^L ® ^L] K 



12; 



(4 + 5i£) 2 + &i£ 7 M L®7 M L 



K. 



22- 



(43c) 



(43d) 



Here Ei[Qj] is the evanescent operator needed as a counterterm to render the 
diagrams C 2 and C3 in Fig. |8] with Qj inserted finite. The colour factors Kij are 



K 12 = K 13 = K l5 
K u = K u = K 16 



1~ 1 
-1 1, 

2 2N ' 

1 N 2 — 2- 

-1 h : 

4 AN 



(44) 



Likewise E 2 [Q 7 ] appears in two-loop diagrams involving physical operators or in 
one-loop matrix elements of E\ [Q 7 ] . 

In our NLO calculation we have kept ai, a 2 , di and hi arbitrary for two 
reasons: First we want to illustrate our findings of P31 in sect. |3.5.1 . Second the 
vanishing of these quantities from physical results provides a non-trivial check of 
our calculation. 

Apart from sect. |3.5.1| we will always state the results corresponding to 

a\ = —8, a 2 = —16, d\ = —8, (45) 
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in order to comply with the standard choice used in [T5|, [If], |2T], [24], ^] . Since 
NLO anomalous dimensions and matching corrections of physical operators do 
not depend on bi, we do not give a numerical value. Likewise we do not need the 
value of the colour factor K22 ■ 

Let us finally look at the operator Q 7 or equivalently at Qs2 appearing in ( pop 
and (P to introduce a different type of evanescent operator: The Dirac and flavour 
structure of Qs2 = (sd)v-A(sd)v-A^- is Fierz self-conjugate in four dimensions. 
Hence Q$2 differs from its Fierz transform (sd)v-A(sd)v-A^- by an evanescent 
operator. In general we must therefore expect the NLO anomalous dimension 
to be different for Q52 and its Fierz transform. Yet the standard definition of 
E\[Qs2\ with Si = —8 ensures Fierz symmetry to hold at the loop level as well 
(see pl| , [33| ). Hence with the choice d\ = —8 it does not matter whether one uses 
Qs2 or its Fierz transform or any linear combination of them in the calculation. 
This is especially gratifying for the hamiltonian (JJ) below the charm threshold, 
if the non-perturbative methods used to obtain the matrix element of Qs2 do not 
distinguish between Qs2 and its Fierz transform. 



3.2.2 Green's Functions from C^^ 2 

The following sections will deal with the determination of the Wilson coefficients 
and renormalization Z-factors present in (j35|) . To do so we need to know the 
Green's function obtained from C^ S<f2 to order Gp. It reads:0 



T exp 



1 J d ^^qQ 



X , 



-i (H c + H l + H ct ) + O (G 3 F ) , (46) 



where 
H c (x) 

if* (x) 
H ct (x) 



C<2 



E/~i /~i iy— 1 iy — l/7-\bare 



X 



1,1 ',3,3 =+>- 



ri2 

\2 fy(t) 9-l/=)bare 



On (x) 



167T : 



2^S2^S2 Qs2 



<JT p 



(47a) 



(47b) 



EEWi E E z u>zj?k^ (x) + zz^r 



x 



i=+-j=l 



TZij (x) 



/~i ry— l^bare 
7 77 ^7 



X 



(47c) 



3 ^The RHS is looking like the Green's function of a |AS|=2 hamiltonian, the notation is a 
little bit sloppy. 
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Here Oij and IZij denote the bilocal structures composed of two | AS| 
reading 



1 operators 



Yj d v T [Q? ( x ) QT (y) - OF ( x ) QY (v) 

_ Q ucMre ^ Qcu^re ^ + gtm.barc ^ g ™,bare ^ 



(48a) 

' r d D y t [ 2 gr ,bare 0*0 Q] u ' haxe (y) - QT ,haxe (*) Qf ,baxe (y) 



Qr ,bare 0*0 QT ,haxe (y) 



for j = 1,2, 



—% 



d D yT (Ql 



,bare ^cc, bare 



Qi 



(x) Qi (y) 



+Q barc (x) (g™> bare _ Q-- bare 



(v). 

for j = 3, ... , 6. 

(48b) 



t /TTo — 1\ j-i j_ /^cc, bare / \ ,otm,bare / \ ,oim,bare / \/-\(X, bare/ \ i i ■ i 

In (|48aj) the term — Qj {x)Qj {y)~Qi \ x )Qj \V) has been omitted, 
because its matrix element does not contribute to the leading power of m^/M^ 
below the charm threshold. Yet this term is necessary to annihilate the mixing 
of the current- current operators into penguin operators in H c . The absence of 
penguins in H c is an effect of the GIM mechanism. In contrast GIM is broken in 
H ct due to the large mass of the top quark. Further the counterterms proportional 
to unphysical operators are not displayed in (|48a ) and ( |48 b| ) . 

Here the LO matrix elements of (|48a|) and (|48b ) with j — 1,2 correspond to 
diagrams of the type shown in Fig. ||, (|48b| ) with j = 3, . . . , 6 to Fig. |14|. 




Figure 14: Diagram Po in the effective five- and four-quark theory. The white 
cross denotes the insertion of a local | AS| = 1 current-current operator, the shaded 
cross the insertion of a local four-quark penguin operator. The penguin op- 
erator only contributes via its up-type-quark foot, i.e. through the couplings 
(sd) v _ A (cc) v±A and {sd) v _ A (uu) v±A . 
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3.3 Matching of the Standard Model Amplitudes to the 
Effective Theory 

Let us now give the initial conditions for the Wilson coefficient functions at the 
scale Htw, & f which the top quark and the W-boson are integrated out. We 
start with the initial conditions needed for the calculation of G ct . It involves the 
coefficients of all |AS| = 1 operators comprised in Q = (Qi, Q 2 , . . ., Qq) t ■ In the 



NLO the initial condition reads pl| : 
/ n \ 



C {ntw) 



a s (/i 



tw 



An 



In 



V>tw 



vv 




1 




+0(a 2 s ). 

To obtain C (ntw) m the LO one simply drops the 0(a s ) term. In ([49] 

11 



( ^ ) 




iff 




(0) 
723 




(0) 
724 
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(0) 
725 




(0) 

V 726 J 





B l 

~m E {xt (ihw)) 

\E{x t (jhw)) 
-^E(x t {ntw)) 
\E(x t (t-kw)) 



B 1 
E(x) 



11 

y 

2, 

— m x + 
3 



B 2 = - 
x 2 (15 



2N 1 

- 16x + Ax 2 ) 



(49) 



(50a) 



6(1 -x) 



Inx + 



»(18-ll*-*»)_2 

12(l-x) 3 3^ ; 



Bi, B 2 and E(x) are scheme dependent, the expressions in ( |50a| ) and ( J50b| ) are 
specific to the NDR scheme and the definition of the evanescent operators given in 



,(o) 



denote the 



(|45j) . The \n(fi t w / Mw) term in (|49D allows for ^^M^, here the 7^ 
elements of the anomalous dimension matrix of the |AS| = 1 operators summarized 
in appendix |C.1| . It is important to note that Qi, Q 2 collectively denote the 
operators Qf, Q 2 l with different flavour quantum numbers k,l = u, c. Note that 
C (ptw) m (HD does not depend on the number / of active flavours, so there is 
no difference whether we match to an effective five-flavour theory or directly to 
an effective four-flavour theory. 

For the forthcoming solution of the RG equations we will also need the diagonal 
basis for the current-current subset of the operator basis, which we have already 
introduced in (EST). The first two rows of (19) then translate into 



tw j 



1 + 



a s (fx 



tw j 



An 



, fkw (o) . D 
In ^^7± + B± 



Mi 



w 



0(a 



with 



B± = B 2 ± Bi, 



(51) 



(52) 
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,(o) 



of the operators Q± can be found in (|154 ). 



and the anomalous dimensions 7± 

We are now in the position to calculate the initial values for the |AS| = 2 
Wilson coefficient C 7 . This is done by comparing the Green's function for the 
|AS| = 2 transition obtained in the full SM with the same quantity obtained in 
the effective five- flavour theory in (|46|). In the SM expression of G ct to order a® 
(cf. ([7a]) and (0)) there is a large logarithm. Therefore the LO matching can 
be done solely with the large logarithm, the NLO matching then requires the 
In part. This can be seen from Table |lj, where one simply has to set n = to 
see the term which is used for the matching in a specific order. It is therefore 
sufficient for both the LO and NLO matching to consider the effective theory only 
to order a®. Since the initial value of the | AS| = 1 coefficients are of order a s 
for k = 1,3,4,5,6 and C k = 1 + 0(a s ) for k = 2,+,- (cf. ([49]) and (§]])), the 
matching at the scale \i = \i t w reads 



2iG ct{0 \fi tw ) + 0(x 2 c \nx c ) 



1 



Gl 



(^ +2 (/w)> (0) + (K- 2 (mw)) w + C 7 {fx tW )(Q 7 {fi tW )} 



<# C W)> 



(o) 



(0) 



,(0) 



(53) 



Here 



has been used. The diagram of Fig. H yields 



while (72._ 2 ) (0) = 0. With fl7a|) and one easily finds 



(0) 



tw, 



16vr 2 



-8 In 



-2 



(QS2) 



(0) 



(54) 



tw, 





47T 



-81n^+4F(x t bH W )) + 2 



Mi 



w 



in LO 
in NLO 



(55) 



where F(x t ) is the top dependent part of S(x c ,x t ) defined in fllPp . The factor a s 
originates from the special definition of Q 7 in (|33|) . Note how the large logarithm 
In x c in ([5c]) is split between the Wilson coefficient C 7 and the matrix element in 
^3|). Again the NLO result in (]54f) and (]55|) is specific to the NDR scheme with 



Next we discuss the other two flavour structures described by G and G c . For 
G l the situation is quite different. Here the SM amplitude immediately has to 
be matched to an effective theory containing only Qs2- The terms needed for 
this matching can again be read off from Table [I] if one sets n = 0. In the LO 
the 0(a®) term of the SM amplitude and of the effective theory matrix element 
is sufficient, i.e. Fig. | and Fig. [Q while in NLO one needs the 0(a]) parts, i.e. 
Fig. |3| and Fig. f|. One then easily extracts the initial condition of the Wilson 



coefficient Cg\ in NLO as 



M w 



S (x t (ntw)) + 



tw j 



4:71 



k l (x t (fitw) , mw) 



(56) 
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where S and k* have been defined in (|9a|) and ( |20c| ) respectively. The LO expres- 
sion is simply obtained by dropping the a s term. Note that the term in square 
brackets in (|56| ) precisely equals the one in (|17a| ). 

The simplest case is G c . Since there is no large logarithm in the SM amplitude 
to order a° due to GIM suppression, we expect the Wilson coefficient of Q-j to 
vanish. We can check this statement explicitly by performing the matching [ |16|j . 
From Table [I] one again reads off the terms required for this by setting n = 0. 
They are the a° terms for LO and the ajln terms for NLO. One therefore in 
LO has to calculate the finite parts of the diagram in Fig. [| with both insertions 
being Q2^s. In NLO one has to do the same with the diagrams in Fig. |15l One 



immediately obtains the result that the double insertions fully account for the 
SM amplitude, there is no room left for Cj 7^ in G c . 







Figure 15: The classes of diagrams yielding the 0(a s 
+, - (|8D and H^; i = +, -; j 



contribution to O^; i,j 



1,2, fl48b| ) in the effective five- and four- 



flavour theory, the remaining diagrams are obtained by mirror reflections. The 
crosses denote insertions of |AS| = 1 current- current operators. Additional QCD 
counterterms have to be included. The result for the divergent parts of D\, i = 
1, . . . , 7 are summarized in appendix ^ ( ( \141a\ ) and Table [|). Diagram D 8 equals 
for zero external momenta. 



3.4 Renormalization Group for Double Operator Inser- 
tions 

So far we have determined the Wilson coefficients taking part in the game at the 
renormalization scale \i t w = 0(Mw,m t ). To obtain them at a scale /i < n t w we 
need to know the RG equation, which governs this evolution. The new feature in 
the calculation presented here is the RG equation for Green's functions with two 
operator insertions. 

The QCD beta function f3 (g), the anomalous mass dimension 7 m and the 
anomalous dimensions for the wave function 7^, are summarized in appendix 0. 
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3.4.1 RG for Single Insertions: A Short Review 

Let us first shortly review the case of single insertions adopting the notation 
of p4| . From /i^£ |As ^ = in (|3l]) one obtains the RG equation 



6 

£ 

3=1 L 



A d 



C< =0 



(57) 



for the Wilson coefficient functions Cj, where 



6 _ d 



(58) 



k=i 



is the anomalous dimension matrix of the | AS| = 1 operators Q/.. The solution of 
floTD is given by 



(59) 



k=l 



with the evolution matrix 

U (/i,/i ) 



T g exp 



(60) 



Here T g means that the matrices j(g'), l{g"), • • • in the expanded exponential 
are ordered such that the couplings increase (decrease) from right to left for 

g(no) < g(ji) (gM > g(jt))- 

We expand the renormalization matrix Z^ 1 as 



z- 1 



Arc 



a. 



47T 



= Y,-r Z r Un) +0{e). 



(61) 



r=0 



To deal with the evanescent operators Z 1 contains a finite renormalization term. 
From ( |5T] ) the coefficients of the perturbative expansion of 



7 



47T 



Art J 



l w + ■■■, 



are obtained as 



7 V 
7 



2Z 1 " 1 ' (1) + 2eZ ^ 



(62) 



(63a) 



(i) 



4Zr 1>(2) + 2 {Z - 1,(1) , zr 1,(1) } + 2[3oZo m + 0(e). (63b) 
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The LO approximation of fl6"0| ) reads 



exp 



d ■ In 



«s (/fo) 
a s (p) 



with 



7 



(0) J 



2Po 



(64) 



Here a s (/i) is the running QCD coupling constant defined in appendix [C], ( |147| ). 
The NLO expression of ( |60D can be written as 



U (/i,/i ) 



47T 



47T 



(65) 



where J is a solution of the matrix equation (see e.g 



1) 



J + [d, J] 



+ 



2/3o A) 



(66) 



We remark here that it is not necessary to diagonalize (|66T ) in order to solve for J. 
flSE|) simply represents a set of 36 linear equations for the 36 elements of J, which 
are therefore rational numbers, d = dffl and J = depend on the number / of 
active flavours through /3q, 7^°- ) and 7^, which can be found in appendix 0. 

For an operator like Q + , Q_, Qj or Q$2 which does not mix with other oper- 
ators, the matrices in fl64| - |6"6"D reduce to numbers. In the following we will need 



\ at(ji) ) 



d T 7 
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~(1) ~(0) a 

777 777 Pi (Rf7 \ 

-% + -w (67a) 

l£> , 7f ft 
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2A, 2/5 2 



-. (67b) 



The anomalous dimensions of Q±, Q 7 and Qs 2 are summarized in appendices |C.l 
and|Tl2l 



During the evolution from \i t w down to fi c we have to pass the scale at 
which we integrate out the bottom quark. Since the penguin operators Qz,...,q 
explicitly depend on the number of active flavours, we have to take into account 
a matching correction r = p4|: 



Q) M 
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(o) . a, 



+ 



4tt 



Q) (1) (/,) + 0(a s 2 ) 



l + ^r(M) + 0(fl«) 



(0) 



(68) 



The matching for the Wilson coefficients C from the effective five-flavour theory 
to the effective four-quark theory therefore reads 



1 + 



a s (/it 
4rr 



] -5r T (ji b ) + O (a£) 



c [5] M 



(69) 
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with 



5r(fi b ) = r^(fi b )-r^(fi b ). 



[-1] 



(70) 



For our |AS| = 1 operator basis the matrix r = r (nib) can be found in ap- 
pendix |C.1| , (|153c|) . The required r (fib) then reads 



(fib) — r ' /J — l n —7 



.[/] 



_ l n i^(0),[/] 



m 6 



(71) 



Since 7^°- ) and r (mj) are independent of / for the current-current subspace, there 
is no matching correction for Qi, Q2- 



3.4.2 An Inhomogeneous RG Equation 

The local operator counterterms proportional to Z^ij (fi) in C^ S]j=s of Q35|) do not 
influence the RG evolution of the coefficients Ci, but they modify the running of 
C-j. We will discuss this in the following. From ft-j^C^ = in (|3~5|) one finds 
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This can be compactly rewritten as 
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Its perturbative coefficients analogous to (|63"D are found as 
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In (|T5]) we have also included finite renormalization constants with subscript 0. 
Such finite renormalizations appear in general when counterterms proportional to 
evanescent operators must be included such as in our calculation. For a detailed 
discussion see |pl|,|33|. The extra terms in (|75b| ) involving the finite renormal- 
ization constants can be simply included into the calculation by multiplying all 
one-loop diagrams containing a finite counterterm by a factor of 1/2.0 

Here and in the following section |3.4.3| we will present two different ways to 
solve the inhomogeneous RG equation (|73|) . With standard methods to solve 
coupled differential equations one obtains 
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Here U$ and J77 are the RG quantities related to single insertions of Q-j defined 
in 



The QCD coupling constant at scale fx has been labeled g and here 
the arguments of the evolution matrices are not the scales but the corresponding 
couplings. 

The first term in (|76|) is solely related to matrix elements with single insertions 
of Q-j. There are no factors involving J77 here, because the initial coefficient C7 (go) 
starts at order g 2 . 

(|76| ) nicely reveals the structure of double insertions: First the two Wilson 
coefficient functions C t > and C s > independently run down from the initial scale 

to the intermediate scale // with g(^i') = g' . Then they are linked by the 
anomalous dimension tensor to the single insertion coefficient C7, which then 
runs further down to the final scale /i. The integral then performs a summation 
over all intermediate scales fi'. If one wants to solve the integral in (p^), one must 
diagonalize at least one of the two |AS| = 1 evolution matrices yielding quite 
cumbersome expressions. 



4 ) Throughout this paper we implement the MS scheme by absorbing 7E — ln(47r) into the 
measure of the loop integrals. Hence this trivial finite part of the counterterms never appears 
explicitly in any formula. 
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3.4.3 A Compact Mixing Matrix 



For formal analyses like those in ]33| the form of (|76|) is well suited. In a practical 



calculation, however, this solution of the inhomogeneous RG equation is difficult 



to implement. Here we present a simpler way to solve (73). 

The key to observation is that in the double insertion diagrams at least one 
of the two operators always stems from the current-current subspace of the full 
| AS | = 1 operator basis. For this subspace we switch to (Q+, Q-), which has the 
advantage that the Wilson coefficient functions C + and C_ do not mix with each 
other as long as we preserve the Fierz-symmetry of Q± during the renormalization 
process. This is the case for our choice of evanescent operators in (|43a|) [2T|, 



The problem then splits into two independent inhomogeneous RG equations 

V-fctb*) = ^M+mAW^W. (77) 
d/i 

Here the decomposition of C-j(iitw) into C^(fi t w) is completely arbitrary provided 
one satisfies 

C 7 (fJLtw) = (lJL tW ) + (5f (jJLtw) . (78) 

This decomposition is then automatically preserved at any renormalization scale. 
We may now cast the inhomogeneous RG equation ([77]) together with the RG 
equations of the |AS| = 1 coefficients into two 7x7 matrix equations: 




7 +7±! 

~T 
7±,7 




(79) 



where 7 is the 6x6 |AS| = 1 anomalous dimension matrix and C = (Ci, . . . , Cq) t . 
Further j± 7 comprises elements of the anomalous dimension tensor defined in 
© and (|7J): 

7±,7 = (7±lj ; 7±2,7,7±3,7,7±4,7,7±5,7,7±6,7) • (80) 

( |79|) and its solution essentially represent the method used by Gilman and Wise 
in their LO analysis [[K|]. Yet they have used an inconvenient operator basis, 



which contains an operator being linearly dependent on the others. The authors 



of |10j therefore involve 8x8 matrices with a double eigenvalue rather than 7x7 
matrices as in (|79|). Further their bilocal structures are defined differently, so 
that they had to solve four RG matrix equations, while we only encounter two of 
them (corresponding to "+" and "— " in (|79|)). 

Yet in ( |79D these two equations still encode a lot of redundant information, 
both evolutions contain the full 6x6 |AS| = 1 evolution matrix. We can do even 
better and collapse them into a single 8x8 RG equation: 

,±D = f-D (81) 
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with 



7 



/ T 



o o \ 

7+, 7 777 -7+ 

777 - 7- / 



V 7-,7 



\ 



C+ (//) /C+ (jm) 
V Cr fa) /C. (/*) y 



.(82) 



Here //^ (p^/C±j appearing on the LHS of directly evaluates to ([T7|) and 
the RG equation for C± involving y±. Both ( [79|) and (|8~ID can be solved by the 
standard techniques already introduced for single insertions in sect. | 3.4.1 
and (1651). 



sec 



Since the Wilson coefficient D contains the Wilson coefficients of the |AS| = 1 
operators, it receives a matching correction when passing from the effective five 
to the effective four quark theory analogously to (|69| ): 

a s (Hb) 



1 + 



47T 



-5r T (fi b ) + O (pi) 



D [5] (/n) 



(83) 



where the 8x8 matrix 5r is defined as 



Sr(n b ) 



( Sr fa) \ 


v j 



»4) 



and 5r denotes the matching correction in the |AS| = 1 Wilson coefficients intro- 
duced in (|70T). 



3.5 The |AS| = 2 NLO Anomalous Dimension Tensor 



In order to calculate the solution of the RG equation (|8T|) we need to know the 
value of the anomalous dimension tensor y±jj, which governs the mixing from 
double insertions to C7. This tensor is determined from the renormalization factor 
^ii,?) see QT§). 

Z^ 7 is determined from the finiteness of the Green's function —i (H ct ) in ( f47c|) . 
Inserting all the required renormalization factors including the wave-function 
renormalization constant we find for the 0(a°), 0{a\) term of ZHf- 



Z-W) 



ft)" *• 



(0),bare 



(85a) 



Q-, 



(0) div. 



E E 7 l 2 4?iM 



%'=+,- j'=i 



. (0),bare 



27 (1) — 



1,(1) 



ij,7 



(0) 



Zi 



-1,(1) 



ij,7 



~ \ (l),barc 

Qr) 

(85b) 
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where we have used the notation (|6T| ) and (|]) for the expanded Z-factors and 
matrix elements. In flSoD the symbol d = means that only the divergent parts of 
the LHS and RHS need to be equal. 



Now in the LO 



ij,7 



is simply obtained from the 1/e pieces Z 1 



1,(1) 



of (7£y) . These terms are calculated by the evaluation of the diagrams in Fig. § 
and Fig. [14| for i = +,— and j = 1, 2 and j = 3, . . . , 6 respectively. 

For the NLO one has to know the 0(a s ) corrections to (Q7), which are related 

to the matrix elements of Q$2 through the definition ([34]): 



Q- 



(l),bare 



/ (~\\ \ 1 / ~ \ (0) 1T1 2 1 / ~ \ (l),bare 

(2Z« + /3 ) - (g 52 ) ( + (Q S2 ) ( • (86) 



Now the divergent parts of the two-loop diagrams in Fig. and Fig. |l^ including 
the corresponding subloop counterterm diagrams yield the terms in the first two 
lines of ( |85b| ) and the last term in ([86]). 
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P 5 Pe P? P 8 




Figure 16: The classes of diagrams yielding the 0(a s ) contribution to IZij, 
i = +, —, j = 3, . . . ,6 in ( \48b\) in the effective five- and four-flavour theory, 
the remaining diagrams are obtained by mirror reflections. The white crosses 
denote insertions of | AS| = 1 current-current operators, the shaded ones insertion 
of | AS | = 1 four-quark penguin operators. Additional QCD counterterms have to 
be included. The divergences of the diagrams P\, i = 1, . . . , 15 are summarized 
in appendix [3| ( (\141b\ ) and Table The diagrams P\ belong to two different 
species: In Pi - Piq the penguin operator contributes through its up-type quark 
foot (as in the LO diagram Pq of Fig. while in P n - P 15 its down-type quark 
foot is involved — a new feature of the NLO calculation. Diagrams vanishing 
identically are not displayed. 
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The remaining divergences therefore correspond to 



-1.(2) 



U,7 



+ A> + 2Z$) z\ 



-1,(1) 



-1,(2) 



(87) 



For clarity counterterms proportional to unphysical operators have been omitted 



in 



Hence we can simply read off Z 1 



-1,(2) 



from the 1/e divergences of the 

diagrams of Fig. [15] and Fig. [16] after the inclusion of subloop counterterms. 

The diagrams Pi in Fig. 1^ appear in two different species: For i — 1, . . . , 10 the 
penguin operator contributes via its up-type quark foot, i.e. through the couplings 
(sd)v-A(cc)v±A and (scI)v-a(uu)v±a as in the LO diagram Po of Fig. [14J. In 
contrast in Pn, . . . , Pi 5 the down-type quark foot of the penguin contributes, 
i.e. the couplings (sd) V -A(ss)v±A and (sd)v-A(dd)v±A- Naively one would not 
expect Pn, . . . , P15 to be proportional to m 2 , since the subdiagram involving m c 
is proportional to 



sin 



Ld ■ (k^k" 



* V) / ( ^ 



i.e. transverse with respect to the virtual gluon momentum k. If one expanded 
f around D = 4, one would only find a logarithmic dependence on m 2 . 

Yet the second loop integration over k is quadratically divergent yielding a result 
proportional to m 2 . Nevertheless only P12 with insertions of Q3 or Q4 has a non- 
vanishing divergent part. The other diagrams are finite. Diagrams with insertions 
of Q5 or Q 6 vanish altogether. Of course there is no divergence proportional to 
m 2 /e 2 in P^, because the one-loop counterterm diagrams vanish. The finiteness 
of Pi 3 — P15 is related to current conservation and ensures that these diagrams do 
not contribute to the NLO calculation of K — > -rvv performed in [53]. We had to 



include the 1PR diagrams P14 and P15 into the consideration, because the result 
of their 1PI subdiagram is proportional to QqEOM defined in (|3D|), but Q q EOM has 
been dropped from the operator basis. 

These unexpected contributions from transverse subdiagrams (|5?f ) has another 
interesting consequence: Dimension-8 operators such as 



Q 9 = sj^LD„d-sD u YLd 

mix into Q7 and even into / m 2 ■ Q7 via two-loop diagrams containing a gluon 
self-energy subdiagram with a c- or 6-quark loop. Hence in an effective field the- 
ory, where gluons can appear in quadratically divergent diagrams heavy degrees 
of freedom (here: the b- and c-quark) of the QCD-lagrangian do not decouple 
anymore. This distinguishes an effective theory with non-renormalizable interac- 
tions (here the four-fermion interactions) from renormalizable theories, in which 
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the Appelquist-Carrazone theorem holds ||35|| . Yet in our case fortunately the 



GIM mechanism ensures that the bilocal structures and IZij of (|48|) do not 
mix into Qg and other physical dimension-8 operators apart from Q7. Hence our 
|AS| = 2 basis is complete. 

Using ( |T5D we then obtain the elements of the anomalous dimension tensor 



-16 
-8 
-32 
-16 
32 
16 



/ 



~(o) 
T-7 



-212 N 
-28 
-456 
-88 



V 



1064 
3 

832 
3 



7^ 



8> 


16 


■16 



276 ^ 
-92 
520 
-216 

1288 



f89a) 



/ 



V 







(89b) 



where we have set N = 3 for brevity. As usual the NLO anomalous dimension 
tensor depends on the renormalization scheme. The result (|89b| ) corresponds to 
the NDR scheme with the definition of the evanescent operators corresponding 
to®. 

In the dia gram P\2 involving the down-type foot of the penguin operator 
contributes —32 to 7±3 j7 and 7±i 7. The results of the individual diagrams can be 
found in appendix Appendix [B] contains j±j for arbitrary N. 



3.5.1 Evanescent Scheme Dependence 

In this section we illustrate some findings of [ 3~3| . In [ 33 1 the transformation rule 
between two anomalous dimension tensors calculated with two different definitions 
of the evanescent operators in ( |4"3"D has been derived. 

In our two-loop calculation we have kept ax, 02, Si and b\ in arbitrary 
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yielding 
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+ — d\ 

4- 32 n 



CL\ ~\~%Ci2 — ^T^l 



3 V 



1$, 



3 -r 3 ui 
-12 +6a! 



1496 
3 



"| a l 



-120 +8ai 

1160 16 



+4ai 

16- 

+4ai 



V" 



3 cix+|a2 + y£Ji 
-128 -4ai-4a 2 -4ai 



.(90) 



Let us first look at the dependence of 7^*7 on a\ and a 2 parameterizing the 
I AS I = 1 evanescent operators in fl43a| - |43b|) : In our case the corresponding formula 
(cf. Eq. (50) of PJ) reads 



7g t7 (a) 



~(0) 

7fc 



•! 7 +Ep-(« , -«)]^7S,7-(9i) 



Here D is a 6x6 diagonal matrix with Di„ 



(i) 



It involves the 



evanescent part — a s /(47r) • Z\ 



-Ei [Qj of the one-loop counterterm to Qi 



needed to render the one- loop diagrams of Fig. [| finite. The normalization in 
(|43|) is chosen such that D is the unit matrix, a and a' are 6x6 matrices defining 
the evanescent operators E\ [Qi] (cf. Eq. (5) of [3"5||). a is easily obtained in terms 
of cti and a 2 from the colour factors in (f43|) and 
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(92) 



where we have chosen iV = 3 for simplicity. For the first term in 
need the current-current part of a in the basis (Q+,Q-)~- 

7 



ID we further 



12 



at, 



a 1 

T' 



a 1 

a 1 
IF' 



(93) 



Here we have tacitly corrected an error in the example at the end of sect. 4 of |33|. 
Inserting (|92|) and (^) into ( |9~ll) with a' = correctly reproduces the dependence 
of 7^7 on ai and a 2 in (|9Up found by our explicit two-loop calculation. 



42 



The dependence on Si is more interesting to study, because it reveals some 
of the subtleties of the Fierz transformation in dimensional regularization: In 
addition to Q 7 we need its Fierz transform 

,2 



g2^2e 



(94) 



with i, j being colour indices. Qj is the Fierz transform of Qs for D = 4, hence 
their difference is evanescent. The one- loop counterterm diagrams involved in the 
calculation of (|90| ) have been accounted for diagram-by-diagram. This effectively 
corresponds to keeping both Q 7 and Q 8 in the operator basis and prevents the 
incorrect use of Fierz symmetry in D- dimensional expressions. For the scheme 
transformation formula we therefore need the one-loop renormalization constants 
in the basis (Q7,Qs)'- 



Z'x 



•1,(1) 



•1,(1) 



•1,(1) 



nF 

+3,7. 
nF 

+j,8 
nF 



-3,7. 
l.(l)l nF 



(0,-2,-16,-4,16,4), 
(-8,-2,0,-4,0,4) 
(0,-2,8,4,-8,-4) 
(4,2,0,-4,0,4). 



(95a) 
(95b) 
(95c) 
(95d) 



Here "nF" means that no Fierz symmetry is used. In the final step to calculate 
the anomalous dimension tensor of the physical operators the |AS|=2 operator 



basis is transformed from (Qr,Q&) to (Q 7 ,Qs — Qt)- Thereby Z 1 



Z[ 



•i,(i) 



nF 



+3,< 



nF 
+3,7 



and 



simply add to Z x 

-,(!) 



1,(1) 



+3,7 



Hence for 7± 7 - 7 the transformation rule |33], Eq. (51)] reads 
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+0,1 

('''I -9-1,(1)1 
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Km" Aim] -(96) 



The extra sum over m compared to [^, Eq. (51)] adds the contributions propor- 
tional to {Qr)^ and (Qs)^ of the diagrams from which 7^-7 is calculated. This 
corresponds to the transformation (Q 7 , Q% — Q 7 ) — > (Q7, Qs) described at the end 
of the previous paragraph. 

The evanescent operator E\ [Qs] is obtained from E x [Q 7 ] in (|43c|) by replacing 
K12 with K\\. The remaining ingredients of ( |96| ) are 



1,(1)1 
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-1,(1)1 
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(97a) 
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(97c) 
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(98) 
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3 

7 -A) 

Inserting (p5|), (|97|-|9"8|) and ( |153a|) into (|96|) with a' = indeed reproduces the 
correct dependence on oi. Note that for j > 3 only the first term in (|96|) is 
nonzero. From e.g. 7147 it is easy to see that one must distinguish Q 7 and Q s to 
derive the correct result. 

The Wilson coefficient Cj depends on Si, too. We find for its initial coefficient: 



An 



-81n^ + AF(x t (fi t w))-(Q + ~ai] 



M. 
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(99) 



Here the dependence on Si enters through 



(^ +2 (^)) (0) 



HW , 



167T 2 

which coincides with fl5"4|) for Si = 
general transformation law: 

Cr(ntWj a i) — Qi^tw i a i) 

E E E [*r w 

fe=+,- j=l l,m=7 



8 In {m c /fit W ) + (6 + Si)] (g S2 ) {0) , (100) 



-8. With the methods of [331 one derives the 



Ck (fJ'tw)Cj (mw) [a' lm — ai r 



(101) 



With the leading order Wilson coefficients Cj^\fi t w) from ([49]) and (|51"|) one easily 
reproduces ( p9|) from (|101[) . 

Finally we discuss the dependence on b\. In [|33| it has been proven that the 
NLO anomalous dimension tensor does not depend on b\. We have verified this 
for our result ( |9"0"D . Nevertheless the individual contributions in ( [75 b| ) to the first 
two components of 7^7 depend on b\. As remarked in |3^] the additive terms 



to 4 



1,(2) 



fcn,7 



in ( |75b|) are automatically taken into account, if one inserts the 



counterterms proportional to evanescent operators with a factor of 1/2 rather 



than 1. In the actual calcu 
A = 1 we have obtained 4 



ation we have inserted them with a factor of A. For 

1,(2)1 



which has been found to depend on b\. 



Setting A = 1/2 in our result yielding 7±j7 makes the coefficient of b\ vanish 
(cf. appendices [S] and |B|). 
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3.6 Should One Sum lnm t /M w ? 



In the calculation of the initial conditions of the Wilson coefficients in sect. ^3 
the top quark and the W-boson have been integrated out simultaneously. This 
procedure is sometimes criticized, because it neglects the RG evolution between 
the scales /i = m t and \i = M w . In |3f| this evolution has been investigated 
in a LO analysis, in which the top quark and the W-boson are integrated out 
separately. While no effect has been found for 7/3, the correction to 772 is claimed 
to be of the same size as the NLO correction calculated in | I5| . Let us therefore 
look at G l defined in (f|) in some detail: The effect of the RG evolution between 
ji = m t and /i = M w is to sum (a s \nx t ) n , n = 0, 1, . . .. The corresponding terms 
for n = and n = 1 are also contained in the LO and NLO results of (|9a| ) and 
( |20c|) . The smallness of both a s (M w ) and \nx t casts doubt on the necessity of 
this extra RG evolution. And, more importantly, any RG summation of In xt is 
accompanied by an OPE: 



G\x u n) = Xf 



p2 



J2C k (m t ,fi)-(O k }(M w ,ri 

. k 



+ c j( m t, /i)C fc (m 4 , /i) 



(S jk ){M w ,fi) 



(102) 



where Sj k are bilocal structures composed of two | AS| = 1 operators describing the 
coupling of an s and d quark to two W-bosons. ( |102|) corresponds to an expansion 



of G l in inverse powers of Xt with higher powers of 1/xt corresponding to increasing 
dimensions of the operators. Hence the price to pay for the summation of the 
small logarithm lna^ to all orders is the inclusion of just a finite number of terms 
in the expansion of G l in l/x t . With the results ( pa| ) and (|20c| ) we can check the 
convergence of this expansion. For x t = 4.6 one finds that even the inclusion of 
the first seven terms in 



_,, N x t 3 . 9 91ns t 15 91nx t 

S{Xt) = J + 2 lnXt -4 + ^2x--4x- t + 



x\ 



4x+ 



In Xt 



.1 / 



(103) 



corresponding to the inclusion of |AS| =2 operators up to dimension 12 in a NLO 
calculation in ( |102j ) results in a 10% error in S(x t ). This error is larger than the 
size of the x t o? s In 2 Xj-term included by the RG summation between m t and M\y. 
Moreover the expansion of k l of ( |20c|) in l/x t shows no convergence at all. Hence 
the RG evolution between /i = m t and /i = Mw is not only unnecessary, but x t 
is simply too small to allow for a meaningful OPE. 

Now in |36[] the summation of lna^ has been tried by the methods of 0, which 
first applies the OPE to the | AS| = 1 substructure of G l and then circumvents the 
calculation of the operator mixing into | AS| = 2 operators by the extraction of the 
relevant logarithms from the |AS| =2 loop diagrams. Yet there are some mistakes 
in the analysis of [3D|: For example the operator basis (Q+, Q 1 ^) has been used for 
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the | AS | = 1 transitions, which is equivalent to shrinking the W- lines to a point as 
in Fig. ||. This is the appropriate method for the case m t -C M w . For m t 3> M w 
one has to shrink the top quark lines in Fig. ^| instead. The expansion of S(x) for 
large x in (|103|) and the one for small x in (|9t]) are obviously different. Further 
the authors of [BO] have not realized that every power of Xt in G l requires different 



operators in (|102| ) with different anomalous dimensions. Hence the results of |36 
are incorrect. 



4 |AS|=2 Transitions at the Charm Threshold 
and Below 

In this section we will eliminate the charm quark as a dynamic degree of freedom 
and describe the physics of the |AS| =2 transition with an effective three-flavour 
lagrangian. The necessary steps are as in the preceding section: 

i) Match the effective four-flavour theory to the effective three-flavour theory 
at the renormalization scale \i c . 

ii) Perform the RG running below /i c . 



4.1 Matching to the Effective Three-Quark Theory 

After integrating out the charm quark all dependence on m c belongs to the Wilson 
coefficients. This implies that the term involving Q7 in (|35| ) has to disappear from 
the effective lagrangian, because Q7 contains m c in its definition (|33|) . Further 
the I AS I = 1 operators are neglected in the new effective lagrangian, because the 
matrix elements of double insertions of these operators are at most proportional 
to m 2 s rather than rr? c . We have already neglected such terms in all preceding 
steps. 

Therefore the new effective lagrangian to describe the physics below fi c reads: 



^ ~ 16^ 



Sbare 



W (a*) + A?Cg (//) + X c X t C [ s f (/*)] Z s l fr) Q\ 2 . 

(104) 

This lagrangian already resembles — H^ as<f2 introduced in (|1]). For the matching 
we have to set the Green's function (|46D and the one derived from (|104|) equal at 
the scale /1 = /i c . 

Let us start the matching procedure with G ct : With ( |47c| ) and (|104|) C 7 (/i c ) 
is calculated from 

E £<3 (jjl c ) Cj (/i c ) (n^) (jjl c ) + c? (/i c ) (q 7 ) (ji c ) 

i=+-j=l 

= ^- 2 C^(^)(Qs2)(fic). (105) 
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CV(/i c ) is already nonzero in the LO due to admixtures from Ci- Recalling the 
inverse power of g 2 in the definition (|33|) of Q 7 one identifies the LO in (|105|) with 

(ct) 



the order a s l . Hence Cg 2 receives a contribution 



2 a s (/i c 



-C 7 (/i c 



in LO. 



(106) 



Since the double insertion diagrams Fig. |5] and Fig. [14] as well as their Wilson 
coefficients are of order a°, they first contribute to order a®, i.e. in the NLO. 
Hence (|106|) is fully sufficient for the LO matching. For the NLO matching we 
define coefficients by 



Then the NLO version of C s c 2 reads 
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(108) 



The coefficients (a*c) i n ( 107 ) are given by the finite parts of the diagrams 
in Fig. [| and Fig. [14]. We find: 
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where the r^'s denote the colour factors 
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(110) 



Note that r^,^ for j = 1,2 depends on the definition of the evanescent operator 
E l \Q 7 ] (cf. (POD ). As usual (|0|) only holds in the NDR scheme. 

If one switches off the RG summation by expanding a s (Htw) an d ® s {Hb) con- 
tained in the LO Wilson coefficient C 7 (Hc) in ( |106| ) around a s (/i c ), one finds 

a s (fi c 



C 7 (fj, c 



Air 



-4m c (hc) In 



9 ■ 
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so that the unfamiliar factor l/a s (/i c ) in ( |106| ) cancels. After inserting ( |1 1 1| ) into 
( |106| ) and the result into ( |104|) the thereby expanded — (C^ S ^~ 2 ) reproduces the 
large logarithm of S(x c ,x t ) in iG ct (cf. (|k|) and (f7a|)). This logarithm appears 
as In (fil/Htw)- K one likewise expands the NLO coefficient in ( |108|) , one finds 
the NLO part F(x t ) of fl9c|) and the two small logarithms In (fif w / M"^) and 
In (ml/nl), which are needed to complete the large logarithm of the LO result to 
lnx c . This shows how the dependence on the matching scales fi c and fi t w cancels 
to the calculated order. 

Let us now shortly discuss the matching for the two remaining cases, i.e. the 
determination of C 



52 



and Cgl in 



104|) . The latter case is particularly simple: 



The corresponding term in the lagrangians (|35|) and (|104| ) is equivalent, the only 



effect for Cgl is t ne transition to the three-quark running of a s . In the case of C 
only the |AS| = 1 operators contribute above /x c , therefore we define coefficients 
dij S2 to parametrize the matching: 



S2 
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Since there is no large logarithm in (Oij (/i)) < -°' ) , the LO matching is performed 
from the finite parts of the diagrams in Fig. [| and Fig. |I|. In the NLO the finite 
parts of Fig. |Tj] and Fig. f| are needed. Expanding dij t s2 in the usual way in a s 
and calculating the required diagrams, we find 



d 



(0) 
ij,S2 



in LO 



(113) 



and 



4&2 M 



T ++ 6 (1 — N) In + 39-33jV-9jV 2 +3jV 3 _|_ 7r 2 -6+6A r +jV 2 -V 3 



12JV 



for 



+ 



T+ 6 (_1 _ TV) In mc( ^ + 13-13iV+3iV 2 -3jV 3 _|_ /R 2 -2+4N-3N 2 +N 3 



127V 



r 6 (—3 — N) ln^ + 



for = (+,-),(-,+) 

-13+3jV+7JV 2 +3jV 3 _|_ /Jl 2 2-6N+5N 2 -N 3 



(114) 



4N 



12N 



for 



in the NLO Here the color factors read 

N + 3 N 
t++ = — - A — , = r_ + = 



N 



(115) 



With (|47aD and flII3ffI3D the NLO Wilson coefficient is found as 



Cs2 (A*c 



i,j=+,- 



d 



(0) 
ij,52 
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a s (/i 



47T 



.(116) 
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4.2 RG in the Effective Three-Quark Theory 



The lagrangian ( |104j ) valid below the renormalization scale /i c only contains the 
single physical operator Qs2- The evolution of the three Wilson coefficients C^, 
j = c,t,ct, for fi < \x c is therefore equal and reads 



a s (p) 



A3] 



1 - J_ 



[3] a s (fj, c ) - a s (p) 



4-7T 



(117) 



where d+ and J+ 1 are the RG quantities for three active flavours defined in (|67b|) . 
The RG evolutions of Qs2 and Q + are equal (see appendix |G.2| , (|155|) ). 
Finally we can express the NLO rji's in (P in terms of the coefficients: 
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c ( sl (aO k (aOF ( i - J- 



Aix 
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Csl Oc) [Ots (f*c)] 
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,[3] 



1 — J. 
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An 



Csf Oc) ["s (Mc)F 



,[3] 



1 - J 



2M^S(x c (fi c ),x t (fi tw )) °< ^ 47r 

The /x-dependence present in ( |117| ) is absorbed into &(/i), which equals 



(118a) 
(118b) 
(118c) 



6 (A*) 



[«- 0)]" 



,[3] 



1+ J. 



[3] gg 0)' 



47T 



(119) 



in the NLO. 

The rji's defined in (|118|) are scheme independent except that they depend on 
the definition of the quark masses. We have adopted the convention of ]n| IEf 
that the running masses in ( |118j ) are defined at the scale at which they are 
integrated out, i.e. m c = m c (/i c ), m t = m t {ntw)- Whenever the r^'s are defined 
such that they multiply S (x c (m c )), S(x t (m t )) and S (x c (m c ) ,x t (m t )) in the 
effective hamiltonian (|l|), we mark them with a star: i]*. For example 



77* S (x c {m c ),x t (m t )) 
We will further use 

m* = m; (mA and 



r] 3 S (x c {n c ),xt (fitw))- 



x A = xArrii 



for i = c,t. 



(120) 



(121) 



The result for in ( |119[ ), however, is scheme dependent through . The 
scheme and scale dependence of 6(/x) must cancel with that in the hadronic matrix 
element of Qs2- 
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5 The Final Result 



In this section we summarize the result and sketch our checks of our NLO calcula- 
tion of 773. Further we give an approximate formula for quick implementations of 
7/3 in phenomenological programs. We close the section with the NLO expressions 
for T]x and 772- 



5.1 The Final Result for 773 in the NLO 

Combining ( 118d.) and ( 108 ) we obtain 



2S (x c (Hc),Xt (Vtw)) 
2vr 



J3] 



1 — J. 



[3] at s {n c 



47T 



/ \ c 7 {n c ) + r ihS2 ivc) Ci (// c ) Cj (/i c ) 

a s KPc) i=+ _ j=1 



(122) 



The Wilson coefficient functions at the renormalization scale fi c , which are needed 
here, are obtained from those at the scale n t w by 



C 7 (/i c 
Cj (/i c 

C± (/i c 



C+ (/i c ) D 7 (// c ) + C_ (// c ) D 8 (/i c ) , 
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(123a) 
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a s (/J, b 



a s (ntw) 
An 



D (ii tw ) , (123d) 



and D(fitw) is obtained from the initial conditions C± (fitw), C (fitw) and 

C 7 (fitw) with the help of (182|) . Further C 7 is split into Cf according to (|78|) . 

The matrices <i = 3-^ and J = encode the 8x8 anomalous dimension matrix 
7 defined in ([82]): 



7 



(o) J 



J 



7 



(i) J 



2/5o A) 



(124) 



50 



in analogy with fl63|) and fl6"6|). 

We emphasize that one should consistently remove terms of order a 2 s in (|123 ) 
and terms of order a s in ( |122|) , because they do not belong to the NLO. 

In Table || we summarize the equations, in which the initial conditions for the 
Wilson coefficients defined at the renormalization scale [i t w as well as the other 



ingredients of (|122|) and (|123| ) can be found. 
Finally rj^ is obtained from r] 3 in ( |122j ) by 



V3 



S ( x c (/i e ) ,x t (ntw)) 
S(x*,x$) 



(125) 



with x* defined in ( |121| ). The exact result for rj^ is independent of fi c , fib and 
fitw- The dependence of our NLO result on these scales will serve as an error 
estimate in sect. Q 



D (fitw) 






C7 (fitw) 


r ij,S2 (/ic) 


S (x c ,x t ) 


d, J 




(ED, © 




GOD, dD 




» (0) 


(124)) 



7 


~(0) ~(1) 
7±,7> 7±,7 


d 






(0) (i) 
7± , 7± 


7 (o) j 7 (i) 


~(0) ~(1) 
777 , 777 


5f 


5r 


r t/l 






( 


67b 


) 


(p4j) 


(153) 


(P7|) 


(H 




(HI) 





A), A 


(fL50aD 





Table 3: The equations in which the ingredients of t \122Q and ( \123{ ) are defined. 



5.2 Analytical Checks 

We have performed several checks of our NLO result in ( 122 - 125]) 



i) The NLO anomalous dimension tensor j±jj (|89b|) has been found indepen- 
dent of the infrared structure of the diagrams in Fig. [15| and Fig. [16], i.e. of 
the small quark masses m s , used as infrared regulators. 

ii) We have kept the gluon gauge parameter £ arbitrary. It has vanished from 
7±j,7 after adding the contributions of the diagrams with their correct com- 
binatorial weight. Further we have checked that £ vanishes from the 1/e 2 



terms in (|87|) after subtracting the ^-dependent term involving Z, 



(i) 



iii) Another check has been provided by the well-known fact that the 1 / e 2 -part 
of the two-loop renormalization constant is related to the one-loop Z-factors 
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or equivalently to the LO anomalous dimensions involved (see e.g. pTJ). We 
have confirmed the corresponding relation for our case: 



-1.(2) 



±jj 



~(o) 

777 



2/?o + 7i 0) 



~(0) 

7±j,7 + 



2^ ljj'l±j',7- 
j'=i 



(126) 



iv) ln(m c / //) / e-terms have disappeared from the sum of two-loop diagrams and 
counterterm diagrams. 

v) The dependences of the final result for 773 on the matching scales /i t wi A*b 
and fi c cancel to order a s . 



vi 



If one expands the final result in powers of a s , one recovers the terms 
proportional to a° In °ln u In 2 I In 1 x c of the result without 

RG improvement in ( [Taj) and ( 17b ). The term proportional to ajln°x c 
cannot be obtained, because it belongs to the NNLO, see Table |TJ The 
expansion of ^3 in terms of a s (/j, c ) reveals how the coefficients of the leading 
and next-to-leading logarithms of k ct in (|20a|) are related to the ingredients 
of the RG calculation: The nonvanishing contributions to the LO term 
proportional to x c In 2 x c are 



16 7m 7+2 < 7 + 32 



lvrv o) +^ (o) 



7)^ + 



16 



7+ 7+2,7 



(127) 



appears here, the remaining part 
Likewise the coefficient of x c In x c 



Note that only the second row of 7^) 
only contributes to higher orders in a s . 
is found to involve the same LO quantities as ( 127 ), the NLO matching 
corrections r ±ji52 (/i c ), F(x t ), B+, E(x t ) (see (|TU7|), (|TDp, (0) and @) and 
the elements 7+2,7 of the NLO anomalous dimension tensor. In addition all 
terms related to the penguin operators sum to zero in the LO and NLO 
part of k ct , e.g. in ( |127| ) one finds 



E 

J=3 



(tS, + tW = o. 



Vll 



The initial condition for C-j in (|55|) as well as the anomalous dimension 
tensor 7± 3 - 7 in ( ^9b| ) depend on the definition of the evanescent operators 
(|43|) . We have checked in sect. |3.5.1| that this dependence is in accordance 
with the theorems of p3[, so that the final result is independent of the 



choice of the evanescent operators. 

We remark that all these checks are not sensitive to the results of diagrams Pn- 
P15, which are gauge independent and have no l/e 2 -divergences. 



52 



5.3 An Approximate Formula for 773 



Since the numerical implementation of 773 in ( 122|) and ( |123|) is quite cumbersome 



for phenomenological studies, we present a simple approximate formula for this 
quantity. 

Such an approximate formula is motivated by the following observations, which 
are derived from the numerical study in sect. || 

i) Variation of fib in the interval \i c < fib < fJ>tw changes the result for 773 on 
the permille level. 

ii) The contribution of the penguin operators Q3, . . . , Qe is of the order of 1%. 
This is so, because penguin effects enter G ct in the order a 2 s rather than a s 
(see vi) in sect. |5\ 



We can therefore simply switch off the penguin operators and further set fib = 
fJ>tw = 0(Mw,fn t ), thereby neglecting any effects from the effective five flavour 
theory. This yields 
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M^r 21 -x t (fi tW ) 
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_ « s (pw) A — A^ 



128) 



1 

12 , ,6 , .24 



A ++ = A%, A + ~ = A~£, A— = A-%. (129) 



We have kept the dependence on the scales \i c and fitw because of their importance 
for our error estimate in sect. |5|. 
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The accuracy of ( |128| ) is 1% with respect to combined variations of /i c , fj, t w, 
m c and m t in reasonable intervals of the parameters. For extremely high values 
of \itw and extremely low values of fi c the precision reduces to something like 2%. 

We can in the same way derive an approximate formula with Hb = [i c = 0(m c ), 
thereby performing the RG evolution in an effective five-flavour theory, but this 
is considerably less accurate than ( |128| ) in the whole parameter space. 



5.4 r]i and 772 



Let us now shortly summarize the results for r/i |15j and rj 2 ||16|| . The former 
is obtained by combining ( |118a| ), (|116|) and (|123c|) . The latter is constructed 
by evolving the Wilson coefficient Cgld^tw) in (EH) down to the scale /i c and 
inserting this into ( |118b| ). One finds 
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where terms not belonging to the LO and NLO have been consistently removed. 
From (|130j) it is immediately clear that the masses enter as m c (/x c ) and m t (fHw)- 
Table ^ refers to the equations in which the various quantities entering fll30|) are 
defined. 

Finally r/* and rfc are defined analogously to r}\ in ( |125| ): 



Vi = Vi 



r} 2 = rj2 



S(xt) 



(131) 
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Table 4: Tie equations in which the ingredients of rjx and r]2 are defined. 

6 Numerical Results 

This section is devoted to the numerical analysis of the results derived in the 
preceding sections. We will present the dependence of n 3 and 773 on its various 
physical parameters and on the renormalization scales at which particles are 
integrated out. Recall the relevant part of the effective low-energy hamiltonian 
from (H): 

RCt M = 7?2- M w 2AA 773 S (x c (X) , x t (ji tw )) b (jj) Q S 2 (ji) 

Ql 

2X c Xt Vt S (z* xt) b (fj.) Q S 2 (ji) ■ (132) 



16vr 

Here 773 and 773 depend on the scales fi t w, and /i c . Further they are functions 
of the masses m t , m c and of Aqcd- To establish a starting point let us pick a 
basic set of input parameters 

m c (m c ) = fi c = 1.3 GeV, A^ = 0.31 GeV (A{g D = 0.15 GeV), 

fit = 4.8 GeV, M w = 80 GeV, (133) 

tnw = 130 GeV, m t (m t ) = 167 GeV. 

In the following Aq C d is always understood to be defined with respect to four 
active flavours, the corresponding quantities in effective three- and five flavour- 
theories are obtained by imposing continuity on the coupling a s at \i c and 

The values of Aqcd quoted in ( |133|) require a comment: The world average 
for a s (M z ) = 0.117 [[37]] corresponds to Aj^g = 310 MeV for /jj = 4.6 GeV and to 



Aj^g = 315 MeV for fif, = 5.0 GeV. The LO Aqq D , however, differs from A^ by 
an overall /i-dependent factor. If one equates the LO coupling and the NLO MS 
coupling constant at the scale /1 = Mz = 91 GeV, one finds that A^g = 310 MeV 
corresponds to Aq^, d = 110 MeV. If the matching relation is imposed at the low 
scale fi = 1.3 GeV, one finds Aqq D = 180 MeV. This shows one contribution to 
the error bar in LO calculations. 



The value for 773 corresponding to the set ( |133| ) reads: 

^ LO = 0.365, r]f LO = 0.467. (134) 



^Threshold corrections appearing for /i q ^ m q are numerically negligible. 
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Hence the NLO calculation has enhanced 7/3 by 27%. From the difference of 0.102 
between the two values in ( [[34] ) 0.022 originates from the change from the LO 
to the NLO running a s . The smallness of this contribution is of course caused 
by the adjustment of Aqq D to fit the NLO running coupling as described in the 
previous paragraph. The explicit 0(a s ) corrections from the NLO mixing and 
matching contribute 0.080. 

Let us list the dominant sources of the enhancement: At the initial scale 
fitw the magnitudes of the coefficients C 2 , C + and C_ are of the order 1, while 
all other coefficients are almost negligible. The RG evolution from fi t w to fi c 
enhances the coefficient C_ by roughly 75% because of the negative sign of the 
anomalous dimension of while the coefficient C + is damped by 25%. Now the 
penguin coefficients C 3 to C 6 are still negligible at /i = /i c , only Ci(/i c ), C 2 (/i c ) and 
CV(/-i c ) are important. In the matching at /i c the contribution of C + (/i c )Ci(/x c ) 
numerically cancels the one of C_( / u c )Ci(/i c ), because the RG damping of the 
former is compensated by a larger colour factor r +1 = 2 having the opposite sign 
of r_i = —1. Finally CV(/z c ) ~ 0.7 has become large due to the RG admixtures 
from C 2 . Hence only C_(/i c )C 2 (/i c ), C + (/i c )C 2 (yU c ) and CV(/z c ) are important. 
In the LO only CV(/ia c ) enters r/ 3 . Keeping only CV(/x c ) in the NLO expression, 
however, overestimates the NLO enhancement by a factor of roughly 1.5, because 
^(/^c) contributes with a negative sign to 773 (for the standard definition of the 
evanescent operators in (|45|)). 

Let us further quantify the influence of the penguin operators Qz,...,q- If one 
neglects them completely, one obtains ^* NLO ' n P = 0.472 with the set in ( |133[ ), i.e. 
their contribution is of the order of 1%. This strong suppression serves as a major 
motivation for the derivation of the approximate formula without penguin effects 



in sect. 5.3 



6.1 Scale Dependence of 773 

Let us now analyze the scale dependences present in 773 in detail. In (|132|) 773 
multiplies S(x c (fi c ), x t (fitw)), which is scale dependent. Ideally their product is 
scale independent. Therefore 773 turns out to be much more useful in the discussion 
of scale dependences, because in ( |132| ) it multiplies S(x*, x*) and other quantities 
which are independent of fi c , fib and fi t w- So 7/3 should essentially behave flat 
with respect to variations of the scales, any remaining dependence may serve as 
a measure of the accuracy of the calculation. 

We now have to discuss the scale dependence of 773 associated with the vari- 
ation of the three scales fi t w, and /i c . The one related to fib appears to be 
extremely mild. This is due to the fact that no diagrams containing internal 
bottom quarks contribute to the |AS| = 2 process in order a s . The only places 
where //& enters are a) the running of a s , b) the NLO matching condition of the 
|AS| = 1 penguin Wilson coefficients ( [[53c] ), and c) the anomalous dimensions of 
the |AS| = 1 penguins. The latter two are strongly suppressed as we have seen 
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already in a preceding paragraph. 

Because the //& dependence of 773 is so weak we could even set [i h = \i c = 0{m c ) 
or fib = fhw — O (raw, rat), as we did for the approximate formula in sect. |5~3[ 
The error introduced by this is of the order of 0.1-0.2%. 

Now let us turn to the more important cases, fi t w and \i c . First consider the 
variation of 7/3 with respect to fitw-, which is displayed in Fig. O. Since at fi t w 




Figure 17: The variation of 7/3 in LO 
and NLO with respect to the scale 
fjLtw, a t which the initial condition is 
dehned. The other input parameters 
are given in (\133j). 



Figure 18: The variation of 773 in LO 
and NLO with respect to the scale fi c . 
The range for the latter is taken un- 
physically large to visualize the break- 
down of perturbation theory. 



the top quark and the W-boson are integrated out simultaneously, it is natural to 
choose the interval Mw < Htw < m t for the analysis. In the LO result for 773 we 
find a sizeable scale dependence of 12%. It is almost totally removed in the NLO, 
where we obtain a variation of less than 3% in this interval. This shows that 
it is very accurate to integrate out the two heavy particles simultaneously. The 
strong improvement in the NLO is due to the smallness of lnxt and vindicates 



our argumentation in sect. p. 6 



The situation is not so nice in the case of the variation of /i c , which is displayed 



in Fig. y| We have intentionally extended the range for /i c to the unphysical low 
value of 0.7 GeV to visualize the breakdown of perturbation theory. Varying \i c 
within the interval 1.1 GeV < fi c < 1.6 GeV yields 



0.33 < 77* LO < 0.40, 0.43 < 77* NLO < 0.50. (135) 
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This corresponds to a reduction of the scale dependence from 20% to 14%. One 
reason for the poor improvement is the fact that the NLO running of the mass is 
stronger than the LO one. 

From ( |135| ) one realizes that the scale variation is not always a good measure 
of the accuracy of the calculation: The central value of r/3 NLO does not lie in the 
range quoted for ?73 L0 . Yet in the NLO one may also judge the contribution of 
the uncalculated 0(a 2 s ) terms by squaring the calculated 0(a s ) corrections. This 
leads to the same interval as quoted in ( |135|) , so that we may consider the given 
interval as a realistic estimate of 773. 

6.2 Dependence of 7/3 on Physical Quantities 

Let us now investigate the dependence of 773 on the physical parameters. From 
the smallness of the coefficient Cj at the initial scale one expects r/3 to be al- 
most independent of m* = mt(mt). This statement is confirmed numerically, see 
Fig. [H5|, allowing to treat r/3 as m r independent in phenomenological analyses. 




180 190 

m t (m t ) [GeV] 



1.2 1.3 1.4 1.5 
m c (m c ) [GeV] 



Figure 19: The dependence of rj^ on 
m t (m t ) in LO and NLO. The other in- 
put parameters are given in ( \133j ). 



Figure 20: The dependence of nt, on 
m r (m r ) in LO and NLO. 



The LO result for 773 depends on m* = m c (m c ) sizeably. Yet this dependence 
is washed out nearly completely if one looks at the NLO 773, see Fig. [21]. 

^-|ash j n (|132|) and physical observables are entered by the product 
n^S(x*,x*). It turns out to be almost independent of m t , but shows a very 
pronounced dependence on m*, due to the dependence of S(x*,x*) on m*. 
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We close this section by a look at the dependence of 773 on Aqcd, which is 
plotted in Fig. [21]. It also turns out to be very moderate. 
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Figure 21: The dependence of 7/3 on 
Aqcd in LO and NLO. Actual val- 
ues for a s (Mz) correspond to Aq° d rs 
0.15 GeV and A m « 0.31 GeV. The 
other input parameters are listed in 



Figure 22: The impact of the NLO 
calculation of 773 on the unitarity tri- 
angle. The parameters in the plot 
are m t (m ( ) = 165 GeV, \V cb \ = 0.041, 
\V ub /V cb \ = 0.08, m c (m c ) = 1.3 GeV 
and B K = 0.7. While the NLO result 
1]% = 0.47 permits two solutions for 
(p,fj), there is no intersection, if the 
old LO value 7/3 LO = 0.37 is used. 



To summarize the findings of our numerical analysis: 773 behaves essentially 
flat with respect to all physical input parameters and therefore may be treated 
as a constant in all phenomenological applications. The dominant uncertainty of 
our estimate of 773 is due to the scale variation with respect to fi c and we quote 
as our final result: 



*NLO 
V 3 



OAT 



-0.03 
-0.04- 



(136) 



6.3 Results for 77* and 773 

The only coefficient which is sensitive to the physical input parameters is rj*. We 
update the result of \W\ to actual values of ct s {Mz) and m* in Table [|. 



Finally the result for r\\ reads fl5 | 



0.573 



+0.003 
-0.010 



(137) 
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m* 


1.20 


1.25 


1.30 
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a s (M z ) 










0.111 
0.114 
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0.120 
0.123 
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1 0R1 +0.229 
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1 .p. n +0.344 
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1 71 o +0.541 
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Table 5: Numerical results for r]\ in the NLO. The remaining input parameters 
are as in ^13$ ) except for \i m = M w . The errors are estimated by varying /i c in 
the interval corresponding to —0.2 < ln(/i c /m*) < 0.2. The accuracy of the NLO 
result is best in the upper right corner and worst in the lower left corner of the 
table. 



for 80 GeV < jj t w < 180 GeV and 1.1 GeV < fi c < 1.6 GeV. The remaining input 
parameters are taken from ( [L33|) , they affect the result marginally. 



6.4 Impact on the Phenomenology of ek 

Today's key quantity to determine the shape of the unitarity triangle is the well- 
measured value for ex [0,|39|]. With our result for n 3 one can now perform this 
analysis with NLO precision. This has been done first in 0, where also the 
Kl — Ks-mass difference has been investigated. 

Hence we only briefly discuss the impact of our new result for n 3 on the 
unitarity triangle here. A useful tool to parametrize the CKM elements are 
the improved Wolfenstein parameters X,A, ~p, rj, for a precise definition we refer 
to 

The constraint from Ek reads |1 



5.3-10" 



B K A^ (l-p)A^rj*S(x* t )+ V * 3 S(x. 



(138) 



Here Bk has been defined in (0). The three terms on the RHS of ( |138| ) contribute 
roughly 75%, 37% and -12%, i.e. r}\ is most important and r]*, which contains 
the largest uncertainties, is least important. If we look at the impact of the NLO 



corrections, however, we find that the one to 773 affects the RHS of ( 138 ) most, 
because the NLO calculation changes 772 by only -3%, while 773 and 77* are shifted 
by 26% and 77% for a s (M z ) = 0.117. This feature is due to the low scales 
affecting 77* and 773. The NLO correction to 773 has the same effect on the RHS of 
(|138|) as a shift in the non-perturbative parameter Bx from 0.75 to 0.82. Hence 
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the uncertainty of the LO result for 773 had influenced the phenomenology in a 
similar way as the present hadronic uncertainty. 

Finally we mention a special feature of the analysis of ex- (|138|) defines a 
hyperbola in the (p, ?])-plane as a function of Br, m< and the magnitude of the 
CKM element Its intersection points with a circle defined by the fourth 
key parameter \V u b/V c b\ are the allowed solutions for the top corner (p, fj) of the 
unitarity triangle. If these four input parameters are too low, one cannot find a 
solution. Hence ( |138|) encodes a "new physics borderline" in the parameter space 
(see |J): If future determinations pin down the input parameters to too low 
values, the Standard Model will be unable to explain the observed CP violation 
in K° — K°-mixing . Todays central values for \V c b\, \V u b/V c b\,m t and Bk are very 
close to this line. Now the NLO shift to 773 has enlarged the allowed range for the 
input parameters and thereby vindicated the CKM mechanism of CP violation. 
This can be seen in Fig. where the hyperbola ( |138[ ) has been displayed for the 
LO and NLO values of 773. 

7 Conclusions 

We have calculated the QCD short distance coefficient 7/3 of the low energy 
I AS I = 2-hamiltonian in the next-to-leading order (NLO) of renormalization group 
improved perturbation theory. The NLO calculation of the other two coefficients 
rji and 7/2 has been repeated. The three coefficients read 

nl = lMtHl, rf 2 = 0.57US, ig = 0.47^4- (139) 

The coefficients are scheme independent except that they depend on the definition 
of the quark masses in H^^ 2 . The results in ( [L39|) correspond to MS-masses 
m c {m^) and m t {m t ) as indicated by the superscript Only 77^ is sensitive to 
the quark masses and a s , the quoted value corresponds to m c {m c ) = 1.3 GeV and 
a s (Mz) = 0.117 GeV. The NLO values in (|139|) have to be compared with the 
old LO results: 

7/1 LO « 0.74, 7^° ^0.59, 7^ LO ^ 0.37. (140) 

The large differences between ( [1401 ) and (|139|) illustrate the big uncertainty of the 
crude leading log approximation. Conceptually the correct definition of quark 
masses and the consistent use of the QCD scale parameter Aj^g can only be 
addressed in quantities calculated beyond the LO. Especially the dependence of 
773 ■ S(x c ,x t ) on m t enters in the NLO. The dependence of 7/3 on renormalization 
scales and physical quantities has been discussed in detail. 

In the standard analysis of the unitarity triangle ex and the B° — B° -mixing 
parameter Xd play pivotal roles. Since the QCD correction factor of the |AB| = 2- 
hamiltonian is also known in the NLO jUj, this analysis is now completely possible 
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with NLO accuracy ||. The shift in 773 between ( |140| ) and ( |139| ) has a larger 
impact on ex than the NLO corrections to 771 and 772- The use of 773 LO imposes 
an error onto e# which is comparable to the present hadronic uncertainty in Bk- 
The NLO shift to 773 has enlarged the range for (|K&|, \V u b/V c b\,m t , Bk) complying 
with the Standard Model explanation of €k- 

We have further explained the construction of the effective lagrangians needed 
to achieve (|139|) . Here we have described the elimination of unphysical operators 
with emphasis on the discussion of evanescent operators, which impose a new type 
of scheme dependence on the NLO anomalous dimension tensor and the Wilson 
coefficients. We have shown that this scheme dependence complies with the 
results of [33| implying its cancellation in 7/3. Further a compact solution for the 
renormalization group equation for Green's functions with two operator insertions 
has been derived. Finally the results of the two-loop diagrams constituting the 
NLO anomalous dimension tensor have been listed. 
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A Results of the Two-loop Diagrams 



Here we list the results for the two-loop diagrams of Fig. [15| and Fig. ITJJ: 



prs 



(4?T 



4 fc) d\ 



(fc) 



0(e°) 



;i41a) 



k 



7. 



(4tt 



A^ 



k 



= 1 15, £ = L,R. 



(A;) (fc) 

P2L + P 1L + O fo\ 

e 2 e v ' 



, (141b) 



Here the one-loop counterterms have been included diagram-by-diagram. The 
result is then proportional to the Dirac structure •y^L^'y^L = (sd)v-A® {sd) V -A 
apart from terms involving evanescent Dirac structures, which are irrelevant for 
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the NLO calculation and are not displayed in ( |141| ). C k s and E r k are colour factors 
listed in Table |] and Table |7|. They depend on the colour structure of the inserted 
operators (Q r ,Q s )- i — L,R is the chirality of the (qq) v ^ A foot of the inserted 
penguin operator. The factor /i~ 2e in (|141|) stems from the definition (|33| ) of Q7. 

In ( 141 b|) an internal charm quark has been assumed, the corresponding result 
with an up quark is zero. The coefficients d\ k ' , in the MS scheme are listed in 
Tables § and || for an arbitrary gluon gauge parameter £. Note the footnote on 
p. |35|. The tabulated values for d[\pig refer to the NDR scheme and the standard 
definition (f45|) of the evanescent operators, while 61 in ( |43d|) is kept arbitrary. The 
evanescent one-loop counterterms have been inserted with a factor of A. For the 
NLO anomalous dimension tensor 7^7 in ( |89 b| ) we need the result with A = 1/2. 

Hence the independence of 7^7 of 61 can be easily verified from Table |] thereby 



illustrating a general feature proven in ||33|| . By comparing those results for D r k 



and P[[ which correspond to diagrams related by the Fierz transformation (such 
as Di and P71J one finds the same results for A = 1/2. This is only true for 
the choice (|4"5| ) adopted in Table || and Table ^. The Dirac algebra has been 
calculated with the help of the computer package Tracer [O] . 



(yQn Qs) 


(Q+,Qi) 


(Q+,Q2) 


(Q-,Qi) 




(Q-,Q 2 ) 


k 












1 

2 
3 
4 
5 
6 
7 


-l+N 2 


-2+N+N 2 


-1+2N-N 2 




-l+N 


AN 

-l-N+N 2 +N 3 


AN 
-2+N+N 2 


AN 

-1+N+N 2 - 


-AT3 


A 

1-N 


AN 
-l+N 2 


AN 
-l+N 


AN 

-1+2N-N 2 




A 



1-N 


AN 

-l-N+N 2 +N 3 


2N 
-2+N+N 2 


AN 

-1+N+N 2 - 




AN 

-l-N+N 2 +N 3 


AN 
-l+N 2 


AN 

-1+N+N 2 - 


N 3 


A 



-1+7V 


AN 
-l+N 2 


2N 
-2+N+N 2 


AN 
-l+N 2 




AN 
-l+N 2 


AN 
-l+N 


AN 
-l+N 2 




4 




AN 


2N 


AN 




(Q r , Q s ) 


(Qi,Q+) 


(Q 2 ,Q+) 


(Qi,Q-) 




(Q2,Q-) 


k 












6 
7 


N-l 
AN 

N-l 
AN 


N 2 -l 

2N 
N 2 +N-2 
AN 


N-l 
AN 

N-l 
AN 






l-Af 
4 



Table 6: The colour factors C r k s of the diagrams in Fig. [15|. Q r (Q s ) is the left 
(right) operator in the diagrams. For k < 5 one has C k s = C k ' . 
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{.Qri Qs) 


(Q + ,Qs,s) 


(Q+, Qa,6) 


(Q-,Q3,s) 




{Q-i Qa,q) 


k 












1 

2 
3 
4 
5 
6 
7 
8 

9,10 
12 


-1-N+N 2 +N 3 


-2+N+N 2 


-1+N+N 2 - 


N 3 


1-N 


AN 

-1-N+N 2 +N 3 


AN 
-1+N 2 


AN 

-1+N+N 2 - 


N 3 


A 



-1+JV 


AN 

-1+JV 2 


2N 
-2+N+N 2 


AN 
-1+N 2 




AN 
-1+N 


AN 
-1+N 2 


AN 
-1+N 




4 




1-JV 
4 

-1+JV 


AN 
-1+N 2 


2N 
-1+N 


AN 
-1+N 2 




AN 
-1+N 

AN 
-1+N 2 


2N 
-2+N+N 2 

AN 
-2+N+N 2 


AN 
-1+N 

AN 
-1+2N-N 2 




AN 
-1+N 2 
AN 

-1-N+N 2 +N :i 


AN 
-1+N 

2N 
-2+N+N 2 


AN 
-1+2N-N 2 
AN 

— 1+N+N 2 - 


N 3 


4 



1-N 


AN 

-1+N 
AN 


AN 
-1+N 
AN 


AN 

-1+N 
AN 




A 

-1+N 
AN 



Table 7: The colour factors ££ s of the diagrams in Fig. In E T k s the first index 
r refers to the current- current operator. 



k 










1 


( m f + m f) (_ 2 £) 


(mf + m^ 


)(-6-3£ + 24A) 




2 


(mf + mf) (-40 


(mf + m^ 


)l(-X + l)bi 




3 


(mf + mf) (6 + 20 


(mf + m 2 j 


) (-l + 3£ + 24A- 




4 


[mf + mf) (-40 


(mf + m 2 


)(6^-|(A-|)&i 


) 


5 


-12m? + 2 (mf + mfj £ 


22m 2 - 3 


(mf + m^) £ 




6 


(mf + mf) (-20 


(mf + m^ 


)3£ 




7 


[mf + mf) (6 + 20 


-11m? - 


23m^ — (mf + mf) 


(se+KA-D&x) 



Table 8: Divergent parts = d^ (mf,m|, X) of the diagrams D k depicted in 
Fig. [1^ according to (141a). Here k labels the diagram. Two different internal up- 
type quark masses m; and mj are involved. In D 5 and D 7 the mass corresponding 
to the upper quark line is m.j. 
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7 

k 


(*0 


(k) 

PlL 


(k) 

P2R 


(fc) 

PlR 




1 




124 


84 


—12 




2 


-12 + 44 


22-64 


12-44 


-16 




3 


-4f 




4£ 


6 




4 


-4£ 


-24 + 64 + 48A 


12 + 44 


-40 + 4£ 




cr 



in i A c 
12 + 44 


-34 - 64 


IO /It 

-12 - 44 


Zb 




6 


12 + 44 


-58 - 64 + 48A 


-4£ 


-18 + 4£ 


u 


7 


-4£ 


12 -6f 


4£ 


-6 




8 


12 + 44 


22 + 64 


-12-44 


-4 




9,10 


-84 





84 







11,13,14,15 
















12 





24 











Table 9: Divergent parts p i£ 
cording to ( \141b\ ). Here k labels the diagram and 
operator present in the inserted penguin operator. 



(A) of the diagrams P|< depicted in Fig. |I6| ac- 



L, R denotes the projection 



B Anomalous Dimension Tensor and Z-Factors 



Here we sketch the relation between the results ( |141| ) for the two-loop diagrams 



and 7±J 7 and 



. We further list the result for j±j 7 for an arbitrary number 



±i,7 



N of colours. 

The 1/e-terms of the two-loop results combine to 



h rs (X) 



k=l 
15 



£ w k (C r k s + C s k r ) (df\ml 0, A) + 4 fe) (0, ml A)) for s=l,2 



E«*2^ (A) 

fe=i 



for s=3,. . . ,6, 



(142) 



where £ = L for s = 3, 4, £ = R for s = 5, 6 and r = +, — . Further u>fc, v k count 
the possible ways of left-right and up-down reflections. On has W\ = w 2 = w 3 = 
= Vi = v 9 = vio = 1/2, Wi = 1/4, and the other weight factors equal unity. 
Now the 1/e-part of the two- loop Z-factor is obtained by 



■1.(2)1 



rs,7 



-h rS (l). 



(143) 



65 



The NLO anomalous dimension tensor is obtained from ( 142|) by 



7« 7 
Irs," 



- ! ; I- 



(144) 



cf. ([73D and sect. [H 

For completeness we list the LO and NLO results in flS9| ) for an arbitrary 
number N of colours: 



(o) 

r+7 



( -4(iV + 1) \ 
-8 
-8(iV + l) 

-16 
8(JV+1) 
16 j 



( 



(o) 

r-7 



V 



A(N - 1) 


8(N- 1) 


-8(JV-1) 




'+,7 



iV- 1 
N 



N-l 

N 



( 6 - 22iV - lliV 2 \ 
12 - UN 
-2 (6 + 22iV+ lliV 2 ) 

-22iV 
2 (4 + KW + lliV 2 ) 
^ 16 (1 + AN) j 

I 6 + 34iV + lliV 2 \ 
-23N 
2 (-6 + 34iV + lliV 2 ) 

-2 (12 + 23iV) 
-2 (-4 + 22iV + lliV 2 




C RG Quantities and Matching Corrections 

In this appendix we collect various quantities needed for the RG evolution. 
The QCD beta function reads 



Pig) 



9' 



16vr 2 



-A) - Pi 



16vr 2 



where the coefficients in LO and NLO are given by 



A) 



[/]_ UN-2f 



P[ f] = f iV 2 - y Nf - 2C F f. 



(145) 



(146) 
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Using this the NLO running QCD coupling constant a s = g 2 /^ir equals 



47T 



/5] 



In In 



/3 ln 



P 2 o In. 



(147) 



in the MS scheme. The LO version of ( 147 ) is obtained by dropping the term 
involving (3% and replacing by Aq°d- 

Further we need the anomalous dimensions for the masses and the fermion 
fields. We define the renormalization constants and anomalous dimensions as 
follows: 



fTlbure—Zmm , 



7 



(148) 



As usual we expand •y m and 7^ in powers of a s \ 



2 \ 2 



16tt 2 



16tt 2 



+ 



X — m,ip. 



(149) 



We need the following coefficients of ( |149| ): 



(0) 



6CV, 



7« = 7« [/] = C F fa 



97 10 
—N-—f 
3 3 J 



(150a) 
(150b) 



To transform masses between different scales we use the NLO evolution equation 



m(fj,) = m (/j ) 



j a s (/x ) - (jit) 



47T 



(151) 



where 



7, 



(0) 



and 



7 



(i) 



2/?o + A) 



(152) 



The LO version of ( |151| ) is obtained by dropping the 0(a s ) term in the right 
bracket. 
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C.l I 

Here we 
|AS| = 1 



AS| = 1 Mixing and Matching Matrices 

collect the important ingredients of the mixing and matching of the 
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(0) 



operator basis in ( pT|) for iV = 3 |24 
( -2 6 
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3 
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81 
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3 
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(153b) 



(153c) 



V 

Here / denotes the number of active flavours. We stress that and r depend on 
the renormalization scheme. The results in ( |153b|) and ( |153cf ) are specific to the 
NDR scheme with the evanescent operators in ( [43a| - |43b] ) defined with a\ = —8 
and ci2 = —16. ( |153c|) is further specific to the 't Hooft-Feynman gauge and the 
choice /i 2 = p 2 where p is a small external momentum used as an IR regulator. In a 
RG improved Wilson coefficient like (B9T) the dependence on the gauge parameter 



and the IR regulator cancels. Results for arbitrary iV may be found in p5 . 

We further give the anomalous dimension of the current- current subspace of 
the I AS I = 1 operators in the diagonal basis Q± for the scheme specified above: 



(0) 
7± 



±6 



iV + 1 



7i 1} 



7V + 1 
2N 



N 3 3 J J 



(154) 
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C.2 | AS| = 2 Anomalous Dimensions 

Here we collect the important quantities for the mixing and matching of the 
| AS | = 2 operator basis in ( |35|) and ( |104j ): Q$2 has the same Dirac structure as the 
|AS| = 1 operators Q±. Since it is Fierz self-conjugate, its anomalous dimension 
is equal to the one of Q + in ( |154j ) 

752 = 7+- (155) 

This result is specific to the definition of the evanescent operator Ei\Qs2\ with 
Si = —8 (cf. sect. |3.2.1| ). Using ([J3|) the anomalous dimension of Q-j can be 
related to 7+: 

777 G?) = 7+ (<7) + 2 7 m (9) + -P (9). (156) 

9 

We therefore get 

# 7 } = 7 f + 2 7 «-2/3 fc , A; = 0,1,... (157) 
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